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ABSTRACT 

A new X-ray selected and X-ray flux-limited galaxy cluster sample is presented. Based on the ROSAT 
All-Sky Survey the 63 brightest clusters with galactic latitude |6ii| > 20deg and flux /x(0.1 — 2.4keV) > 2 x 
10"^"'^ ergss"^ cm~^ have been compiled. Gravitational masses have been determined utilizing intracluster gas 
density profiles, derived mainly from ROSAT PSPC pointed observations, and gas temperatures, as published 
mainly from ASCA observations, assuming hydrostatic equilibrium. This sample and an extended sample of 
106 galaxy clusters is used to establish the X-ray luminosity-gravitational mass relation. From the complete 
sample the galaxy cluster mass frmction is determined and used to constrain the mean cosmic matter density 
and the amplitude of mass fluctuations. Comparison to Press-Schechter type model mass functions in the 
framework of Cold Dark Matter cosmological models and a Harrison-Zeldovich initial density fluctuation 
spectrum yields the constraints fim = 0.121q'q4 and as = 0.96lgi[2 (90% c.L). Various possible systematic 
uncertainties are quantifled. Adding all identified systematic uncertainties to the statistical uncertainty in 
a worst case fashion results in an upper limit flm < 0.31. For comparison to previous results a relation 
as = 0.4317"'^ '^* is derived. The mass function is integrated to show that the contribution of mass bound 
within virialized cluster regions to the total matter density is small, i.e., Ociuster 
masses larger than 6.4+o;g x lO" h^^ Mq. 
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Subject headings: cosmological parameters — cosmology: observations — galaxies: clusters: general — intergalactic 
medium — large-scale structure of universe — X-rays: galaxies: clusters 



1. Introduction 

The galaxy cluster mass function is the most funda- 
mental statistic of the galaxy cluster population. It is 
determined by the initial conditions of the mass distri- 
bution set in the early universe in a relatively straight- 
forward way, since the evolution of the large-scale matter 
distribution on scales comparable to the size of clusters 
and larger is linear and since the formation of clusters 
is governed by essentially only gravitational processes. 
Choosing a specific cosmological scenario provides defini- 
tive predictions about these initial conditions in a sta- 
tistical sense. The overall process of the gravitational 
growth of the density fluctuations and the development 
of gravitational instabilities leading to cluster formation 
is comparatively easy to understand. It has been well 
described by analytical models (e.g.. Press & Schechter 
1974, Bond et al. 1991, Lacey & Cole 1993, Kitayama 
& Suto 1996, Schuecker et al. 2001a) and simulated in 
numerical gravitational TV-body calculations (e.g., Efs- 
tathiou et al. 1988; Lacey & Cole 1994). Even though the 
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slight deviations from the theoretical prescription found 
in recent simulations (e.g., Governato et al. 1999; Jenkins 
et al. 2001) seem to require further refinement in the the- 
oretical framework, these deviations are too small to be 
significant for the current investigation, as will be shown 
later. For the precision needed here these large simula- 
tions therefore support the validity of the model predic- 
tions. Within this framework the observed cluster mass 
function provides the opportunity to test different cos- 
mological models. The tests are particularly sensitive to 
the amplitude of the cosmic matter density fluctuations 
(at a scale of the order of lOMpc) as well as the nor- 
malized total matter density, (e.g., Henry & Arnaud 
1991; Bahcafl & Cen 1992). 

In addition to its importance in testing cosmological 
models the integral of the mass function yields the in- 
teresting information on the fractional amount of matter 
contained in gravitationally bound large-scale structures. 
Using one of the first attempts to construct a mass func- 
tion over the mass range from giant ellipticals to mas- 
sive clusters by Bahcall & Cen (1993) Fukugita et al. 
(1998) obtain a mass fraction rJcroup = 0.12 ± 0.02 for 
2 X 10^2 < M <2 X 10^^ h^^^ Mq (where the mass frac- 
tion is expressed here in units of the critical density of 
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the universe, pc). This result is already close to the total 
matter density in some of the proposed cosmologic;al sce- 
narios. Therefore a precise observational determination 
of the integral mass function is a very important task for 
astronomy. 

Unfortunately the galaxy cluster mass is not an eas- 
ily and directly observable quantity (except for measure- 
ments of the gravitational lensing effect of clusters which 
may play a large role in the construction of mass func- 
tions in the future) and one has to resort to the observa- 
tion of other cluster parameters from which the cluster 
masses can be deduced. 

X-ray astronomy has provided an ideal tool to first de- 
tect and select massive clusters by measuring their X-ray 
luminosity and to secondly perform mass determinations 
on individual clusters through X-ray imaging and X-ray 

spectroscopy. 

In this paper we report the first rigorous application 
of these two approaches for the construction of the clus- 
ter mass function. Building on the ROSAT All-Sky X- 
ray Survey (RASS) (Triimper 1993; Voges ct al. 1999), 
which has been well studied in the search of the brightest 
galaxy clusters through several survey projects (see refs 
in Sect. 2), we have compiled a new, highly complete sam- 
ple of the X-ray brightest galaxy clusters (HIFLUGCS, 
the Highest X-ray FLUx Galaxy Cluster Sample). 

Thanks to the numerous detailed galaxy cluster obser- 
vations performed with the ROSAT and ASCA (Tanaka 
et al. 1994) satellite observatories and accumulated in the 
archives we can perform a detailed mass determination 
for most of these clusters and obtain a good mass esti- 
mate for the few remaining objects. From these data we 
first establish a (tight) correlation of the measured X-ray 
luminosity and the cluster mass. This relation assures 
that we have essentially sampled the most massive clus- 
ters in the nearby universe, which forms the basis of the 
construction of the cluster mass function. 

Previous local galaxy cluster mass functions have been 
derived by Bahcall & Cen (1993), Biviano et al. (1993), 
Girardi et al. (1998), and Girardi & Giuricin (2000, for 
galaxy groups). Bahcall & Cen (1993) used the galaxy 
richness to relate to cluster masses for optical observa- 
tions and an X-ray temperature-mass relation to con- 
vert the temperature function given by Henry & Arnaud 
(1991) to a mass function. Biviano et al. (1993), Gi- 
rardi et al. (1998), and Girardi & Giuricin (2000) used 
velocity dispersions for optically selected samples to de- 
termine the mass function. Here we use a different ap- 
proach and construct the first mass function for X-ray 
selected galaxy clusters based on the RASS using indi- 
vidually determined cluster masses. The mass function 
of this cluster sample is then used to determine the mass 
fraction in bound objects with masses above a minimum 
mass and to derive tight constraints on cosmological sce- 
narios. 



The paper is organized as follows. In Sect. 2 the sam- 
ple selection is described. The details of the determina- 
tion of the observational quantities are given in Sect. 3. 
The results are presented in Sect. 4 and discussed in 
Sect. 5. The conclusions are summarized in Sect. 6. 

Throughout a Hubble constant Hq = 50 h^o kms~^ Mpc 
hso = 1, Om = 1, a normalized cosmological constant 
f^A = ^/{^Hq) = 0, and a normalized curvature index 
rik = 1 — f^m — f^A = is used if not stated otherwise 
(present day quantities). We note that the determi- 
nation of physical cluster parameters has a negligible 
dependence on and Ha for the small redshift range 
used here, as will be shown later. Therefore it is justified 
to determine the parameters for an Einstein-de Sitter 
model but to discuss the results also in the context of 
other models. 

2. Sample 

The mass function measures the cluster number den- 
sity as a function of mass. Therefore any cluster fulfilling 
the selection criteria and not included in the sample dis- 
torts the result systematically. It is then obvious that for 
the construction of the mass function it is vital to use a 
homogeneously selected and highly complete sample of 
objects, and additionally the selection must be closely 
related to cluster mass. In this work the RASS, where 
one single instrument has surveyed the whole sky, has 
been chosen as the basis for the sample construction. Us- 
ing the X-ray emission from the hot intracluster medium 
for cluster selection minimizes projection effects and the 
tight correlation between X-ray luminosity and gravita- 
tional mass convincingly demonstrates that X-ray cluster 
surveys have the important property of being mass selec- 
tive (Sect. 4.1). 

Several cluster catalogs have already been constructed 
from the RASS with high completeness down to low flux 
limits (see refs below). These we have utilized for the 
selection of candidates. Low thresholds have been set for 
selection in order not to miss any cluster due to measure- 
ment uncertainties. These candidates have been homoge- 
neously reanalyzed, using higher quality ROSAT PSPC 
pointed observations whenever possible (Sect. 3). A flux 
limit well above the limit for candidate selection has then 
been applied to deflne the new fiux-limited sample of the 
brightest clusters in the sky. 

In detail the candidates emerged from the following 
input catalogs. Table 1 lists the selection criteria and 
the number of clusters selected from each of the cata- 
logs, that are contained in the final sample. 

1) The REFLEX (ROSAT-ESO Flux-Limited X-ray) 
galaxy cluster survey (Bohringer et al. 2001b) covers the 
southern hemisphere (declination 6 < -1-2.5 deg; galactic 
latitude |6ii| > 20.0deg) with a flux limit /x,iim(0.1 — 
2.4 keV) = 3.0 X lO-i^ergss-^cm-^. 

2) The NORAS (Northern ROSAT All-Sky) galaxy clus- 
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ter survey (Bohringer et al. 2000) contains clusters show- 
ing extended emission in the RASS in the northern hemi- 
sphere {S > 0.0 deg; > 20.0 deg) with count rates 
Cx(0.1 - 2.4 keV) > 0.06 ctss-^ 

3) NORAS II (.1. Retzlaff ct al., in preparation) is the 
continuation of the NORAS survey project. It includes 
point like sources and aims for a flux limit /x,iim(0.1 — 
2.4keV) = 2.0 x lO'^^ (^j-ggg-i ^m-z. 

4) The BCS (ROSAT Brightest Cluster Sample) (Ebel- 
ing et al. 1998) covers the northern hemisphere {S > 
0.0 deg; > 20.0 deg) with /x,iir„(0.1 - 2.4 keV) = 

4.4 X 10~^^ ergss"^ cm~^ and redshifts z < 0.3. 

5) The RASS 1 Bright Sample of Clusters of Galaxies (De 
Grandi et al. 1999) covers the south galactic cap region in 
the southern hemisphere (S < +2.5 deg; bu < —20.0 deg) 
with an effective flux limit /x,iim(0.5 — 2.0 keV) between 
~3 and 4 x 10^^^ ergs s~^ cm~^. 

6) XBACs (X-ray Brightest Abell-type Clusters of galax- 
ies) (Ebeling et al. 1996) is an all-sky sample of Abell 
(1958)/ACO (AbeU et al. 1989) clusters hmited to high 
galactic latitudes |6ii| > 20.0 deg with nominal AGO red- 
shifts z < 0.2 and X-ray fluxes /x(0.1 - 2.4 keV) > 
5.0 X 10~^^ ergss"^ cm~^. 

7) An all-sky list of Abell/ AGO/ ACO-supplementary 
clusters (H. Bohringer 1999, private communication) 
with count rates Cx(0.5 - 2.0 keV) > 0.6ctss-^ 

8) Early type galaxies with measured RASS count rates 
from a magnitude limited sample of Beuing et al. (1999) 
have been been checked in order not to miss any X-ray 
faint groups. 

9) All clusters from the sample of Lahav et al. (1989) and 
Edge et al. (1990), where clusters had been compiled from 
various X-ray missions, have been checked. 

The main criterion for candidate selection, a flux 
threshold 1.7 x 10~^^ ergss"^ cm~^, has been chosen to 
allow for measurement uncertainties in the input cat- 
alogs. E.g., for REFLEX clusters with 1.5 < /x < 

2.5 x 10~^^ ergss"^ cm~^ the mean statistical flux error 
is less than 8%. With an additional mean systematic 
error of 6 %, caused by underestimation of fluxes due to 
the comparatively low RASS exposure times^, the flux 
threshold 1.7 x 10~^^ ergss"^ cm~^ for candidate selec- 
tion then ensures that no clusters are missed for a final 
flux limit /x,iim = 2.0 x 10"-^^ ergss^^ cm^^. 

Almost none of the fluxes given in the input cata- 
logs have been calculated using a measured X-ray tem- 
perature, but mostly using gas temperatures estimated 
from an Lx~2x relation. In order to be independent 
of this additional uncertainty clusters have also been se- 
lected as candidate's if they exceed a count rate thresh- 

■^This has been measured by comparing the count rates determined 
using pointed observations of clusters in this work to count rates for 
the same clusters determined in REFLEX and NORAS. If count 
rates are compared also for fainter clusters, not relevant for the 
present work, the mean systematic error increases to about 9 % 
(Bohringer et al. 2000). 



old which corresponds to /x(0.1 - 2.4 keV) = 2.0 x 
10^^^ ergs s~^ cm^^ for a typical cluster temperature, 
Igas = 4keV, and redshift, z = 0.05, and for an excep- 
tionally high column density, e.g., in the NORAS case 
riH = 1-6 X lO^icm-2. 

Most of the samples mentioned above excluded the 
area on the sky close to the galactic plane as well as 
the area of the Magellanic Clouds. In order to con- 
struct a highly complete sample from the candidate list 
we applied the following selection criteria that successful 
clusters must fulfil: 

1) redetermmerf flux /x(0. 1-2.4 keV) > 2.0x 10^^^ ergss"^ cm~^, 

2) galactic latitude |&ii| > 20.0 deg, 

3) projected position outside the excluded 324 deg^ area 
of the Magellanic Clouds (see Tab. 2), 

4) projected position outside the excluded 98 deg^ region 
of the Virgo galaxy cluster (see Tab. 2)."^ 

These selection criteria are fulfilled by 63 candidates. 
The advantages of the redetermined fluxes over the fluxes 
from the input catalogs are summarized at the end of 
Sect. 3.1. In Tab. 1 one notes that 98% of all clusters 
in HIFLUGCS have been flagged as candidates in RE- 
FLEX, NORAS, or in the candidate list for NORAS II; 
these surveys are not only all based on the RASS but 
all use the same algorithm for the count rate determina- 
tion, further substantiating the homogeneous candidate 
selection for HIFLUGCS. 

The fraction of available ROSAT PSPG pointed obser- 
vations for clusters included in HIFLUGCS equals 86 %. 
The actually used fraction is slightly reduced to 75 % be- 
cause some clusters appear extended beyond the PSPG 
field of view and therefore RASS data have been used. 
The fraction of clusters with published ASCA temper- 
atures equals 87%. If a lower fltix limit had been cho- 
sen the fraction of available PSPG pointed observations 
and published ASCA temperatures would have been de- 
creased thereby increasing the uncertainties in the de- 
rived cluster parameters. Furthermore this value for the 
flux limit ensures that no corrections, due to low expo- 
sure in the RASS or high galactic hydrogen column den- 
sity, need to be applied for the effective area covered. 
This can be seen by the effective sky coverage in the RE- 
FLEX survey area for a flux limit /x,iim(0.1 — 2.4 keV) = 
2.0 X 10~^^ ergss"^ cm~^ and a minimum of 30 source 
counts, which amounts to 99 %. The clear advantage is 
that the HIFLUGCS catalog can be used in a straight- 
forward manner in statistical analyses, because the effec- 
tive area is the same for all clusters and simply equals 
the covered solid angle on the sky. 

The distribution of clusters included in HIFLUGCS 
projected onto the sky is shown in Fig. 1. The sky 



^The large scale X-ray background of the irregular and very ex- 
tended X-ray emission of the Virgo cluster makes the undiscrimi- 
nating detection/selection of clusters in this area difficult. Candi- 
dates excluded due to this criterium are Virgo, M86, and M49. 
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Table 1 
Selection of candidates 



Catalog 




Cx [ctss-i] 


/x [10-11 erf 


;s s 1 cm ^] 


iVci 


Ref. 


0.1 - 2.4 kcV 


0.5-2.0keV 0.64- 2.36keV 


0.1 - 2.4 keV 


0.5- 2.0keV 




REFLEX 




0.9 


1.7 




33 


1 


NORAS 




0.7 


1.7 




25 


2 


NORAS 11 




0.7 


1.7 




4 


3 


BCS 


1.0 




1.7 




1 


4 


RASS 1 








1.0 





5 


XBACs 


1.0 




1.7 







6 


Abell/ACO 




0.7 









7 


Early type gx 




0.7 









8 


Previous sat^ 
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''All clusters from this catalog have been flagged as candidates. 



Note. — Only one of the criteria, count rate or flux, has to be met for a cluster to be selected as candidate. The catalogs 
are listed in search sequence, therefore iVci gives the number of candidates additionally selected from the current catalog 
and contained in the flnal flux-limited sample. So in the case of NORAS a cluster is selected as candidate if it fulfils 
Cx(0.5 - 2.0keV) > 0.7ctss-i or /x(0.1 - 2.4keV) > 1.7 x 10-" ergs s 1 cm ^ and has not already been selected from 
REFLEX. This candidate is counted under Nqi if it meets the selection criteria for HIFLUGCS. 



References. — (1) Bohringer et al. 2001b. (2) Bohringer et al. 2000. (3) J. Retzlaff et al., in preparation. (4) Ebeling 
et al. 1998. (5) De Grandi et al. 1999. (6) Ebeling et al. 1996. (7) H. Bohringer 1999, private communication. (8) Beuing 
et al. 1999. (9) Lahav et al. 1989; Edge et al. 1990. 



Table 2 

Regions of the sky not sampled in hiflugcs 



Region 


R.A. Range 


Dec. Range 


Area (sr) 


LMC 1 


58 103° 


-63 -+ -77° 


0.0655 


LMC 2 


81 -> 89° 


-58 -63° 


0.0060 


LMC 3 


103 ->■ 108° 


-68 -74° 


0.0030 


SMC 1 


358.5 -> 20° 


-67.5 ->■ -77° 


0.0189 


SMC 2 


356.5 -> 358.5° 


-73 ->■ -77° 


0.0006 


SMC 3 


20 30° 


-67.5 -72° 


0.0047 


Virgo 


182.7 192.7° 


7.4 17.4° 


0.0297 


Milky Way'' 


-> 360° (In) 


-20 -^ 20° (6ii) 


4.2980 



''Galactic coordinates. 



Note. — Excised areas for the Magellanic Clouds are the same 
as in Bohringer et al. (2001b), because REFLEX forms the basic 
input catalog in the southern hemisphere. 
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coverage for the cluster sample equals 26 721.8 deg^ 
(8.13994sr), about two thirds of the sky. The cluster 
names, coordinates and redshifts are listed in Tab. 3. 
Further properties of the cluster sample are discussed in 
Sect. 5.1. 

For later analyses which do not necessarily require a 
complete sample, e.g., correlations between physical pa- 
rameters, 43 clusters (not included in HIFLUGCS) from 
the candidate list have been combined with HIFLUGCS 
to form an 'extended sample' of 106 clusters. 




Fig. 1. — Aitoff projection of the 63 HIFLUGCS galaxy 
clusters in galactic coordinates (filled circles). Addition- 
ally shown are 11 clusters above the flux limit but with 
|6ii| < 20.0deg (open triangles). 



3. Data Reduction and Analysis 

This Section describes the derivation of the basic 
quantities in this work, e.g., count rates, fluxes, luminosi- 
ties, and mass estimates for the galaxy chisters. These 
and other relevant cluster parameters are tabulated along 
with their uncertainties. 

3.1. Flux Determination 

Measuring the count rate of galaxy clusters is an im- 
portant step in constructing a flux-limited cluster sample. 
The count rate determination performed here is based on 
the growth curve analysis method (Bohringer et al. 2000), 
with modifications adapted to the higher photon statis- 
tics available here. The main features of the method as 
well as the modifications are outlined below. 

The instrument used is the ROSAT PSPC (Pfeffer- 
mann et al. 1987), with a low internal background ide- 
ally suited for this study which needs good signal to 



noise of the outer, low surface brightness regions of the 
clusters. Mainly pointed observations from the public 
archive at MPE have been used. If the cluster is extended 
beyond the PSPC field of view making a proper back- 
ground determination difficult or if there is no pointed 
PSPC observation available, RASS data have been used. 
The ROSAT hard energy band (channels 52 - 201 « 
0.5 — 2.0 keV) has been used for all coimt rate measure- 
ments because of the higher background in the soft band. 

Two X-ray cluster centers are determined by find- 
ing the two-dimensional 'center of mass' of the photon 
distribution iteratively for an aperture radius of 3 and 
7.5arcmin around the starting position. The small aper- 
ture yields the center representing the position of the 
cluster's peak emission and therefore probably indicates 
the position where the cluster's potential well is deep- 
est. This center is used for the regional selection, e.g. 
^11 1 > 20.0. The more globally defined center with the 
larger aperture is used for the subsequent analysis tasks 
since for the mass determination it is most important to 
have a good estimate of the slope of the surface bright- 
ness profile in the outer parts of the cluster. 

The background surface brightness is determined in a 
ring outside the cluster emission. To minimize the in- 
fluence of discrete sources the ring is subdivided into 
twelve parts of equal area and a sigma clipping is per- 
formed. To determine the count rate the area around 
the global center is divided into concentric rings. For 
pointed observations 200 rings with a width of 15arcsec 
each are used. Due to the lower photon statistics a width 
of 30arcsec is used for RASS data and the number of 
rings depends on the fleld size extracted (100 — 300 rings 
for field sizes of 2 x 2 deg^ — 8 x 8deg^). Each photon 
is divided by the vignetting and deadtime corrected ex- 
posure time of the skypixel where it has been detected 
and these ratios are shimmed up in each ring yielding the 
ring count rate. From this value the background count 
rate for the respective ring area is subtracted yielding 
a source ring count rate. These individual source ring 
count rates are integrated with increasing radius yield- 
ing the (cumulative) source count rate for a given radius 
(Fig. 2). Obvious contaminating point sources have been 
excluded manually. The cut-out regions have then been 
assigned the average surface brightness of the ring. If 
a cluster has been found to be clearly made up of two 
components, for instance A3395n/s, these components 
have been treated separately. This procedure ensures 
that double clusters are not treated as a single entity for 
which spherical symmetry is assumed. For the same rea- 
son strong substructure has been excluded in the same 
manner as contaminating point sources. In this work the 
aim is to characterize all cluster properties consistently 
and homogeneously. Therefore if strong substructure is 
identified then it is excluded for the flux/luminosity and 
mass determination. 

An outer significance radius of the cluster, rx, is de- 
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termined at the position from where on the Poissonian 
l-cr error rises faster than the source count rate. Usually 
the source count rate settles into a nearly horizontal line 
for radii larger than rx- We have found, however, that 
in some cases the source coimt rate seems to increase 
or decrease roiighly quadratically for radii larger than 
rx indicating a possibly under- or overestimated back- 
ground (Fig. 3). We therefore fitted a parabola of the 
form y = mx^ + 6 to the source count rate for radii larger 
than rx and corrected the measured background. An ex- 
ample for a corrected source count rate profile is shown 
in Fig. 4. 

Figure 5 shows for the extended sample (106 clus- 
ters) that the difference between measured and corrected 
source count rate is generally very small. Nevertheless an 
inspection of each count rate profile has been performed, 
to decide whether the measured or corrected count rate 
is adopted as the final count rate, to avoid artificial cor- 
rections due to large scale variations of the background 
(especially in the large RASS fields). The count rates are 
given in Tab. 3. 




10 20 30 40 50 

A2029 rac)ius_a (arcmin) 



Fig. 2. — Cumulative source count rate as a function 
of radius (solid line) for the cluster A2029 (pointed ob- 
servation). The vertical dashed line indicates the outer 
significance radius, rx- The dashed lines just above and 
below the source count rate indicate the l-cr Poissonian 
error bars. 

The conversion factor for the count rate to flux conver- 
sion depends on the hydrogen column density, nn, on the 
cluster gas temperature, Tgas, on the cluster gas metal- 
licity, on the cluster redshift, z, and on the respective de- 
tector responses for the two different PSPCs used. The 
njj value is taken as the value inferred from 21 cm radio 
measurements for our galaxy at the projected cluster po- 
sition (Dickey & Lockman 1990; included in the EXSAS 
software package, Zimmermann et al. 1998; photoelec- 
tric absorption cross sections are taken from Morrison 
& McCammon 1983). Gas temperatures have been es- 
timated by compiling X-ray temperatures, Tx, from the 
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20 40 60 80 100 

EX00422 radius_a (arcmin) 

Fig. 3. — Cumulative source count rate as a function of 
radius for the cluster EXO0422, shown as an extreme ex- 
ample (RASS data). The parabolic dashed line indicates 
the best fit parabola for count rates larger than rx- 
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20 40 60 SO 100 

EXD0422 radius_a (arcmin) 

Fig. 4. — Corrected cumulative source count rate as a 
function of radius for the cluster EXO0422. The count 
rate correction is less than 5 %. 



6 




0.1 1.0 10.0 

measured ctr [cts/sec] 

Fig. 5. — Comparison of measured and corrected source 
count rates for the extended sample of 106 galaxy clus- 
ters. The solid line indicates equality. 

literature, giving preference to temperatures measured 
with the ASCA satellite. For clusters where no ASCA 
measured temperature has been available, Tx measured 
with previous X-ray satellites have been used. The X-ray 
temperatures and corresponding references are given in 
Tab. 4. For two clusters included in HIFLUGCS no mea- 
sured temperature has been found in the literature and 
the Lx{< 2 h^f^ Mpc) - Tx relation of Markevitch (1998) 
has been used. The relation for non cooling flow cor- 
rected luminosities and cooling flow corrected/emission 
weighted temperatures has been chosen. Since the con- 
version from count rate to flux depends only weakly on 
Tgas in the ROSAT energy band for the relevant tempera- 
ture range a cluster temperature fcTgas = 4 keV has been 
assumed in a first step to determine Lx(< 2 h^^^ Mpc) for 
the clusters where no gas temperature has been found in 
the literature. With this luminosity the gas temperature 
has been estimated. The mctallicity is set to 0.35 times 
the solar value for all clusters (e.g., Arnaud et al. 1992). 
The redshifts have been compiled from the literature and 
are given in Tab. 3 together with the corresponding ref- 
erences. With these quantities and the count rates given 
in Tab. 3 fluxes in the observer rest frame energy range 
0.1 — 2.4 keV have been calculated applying a modern 
version of a Raymond-Smith spectral code (Raymond & 
Smith 1977). The results are hsted in Tab. 3. The flux 
calculation has also been checked using XSPEC (Arnaud 
1996) by folding the model spectrum created with the 
parameters given above with the detector response and 
adjusting the normalization to reproduce the observed 
count rate. It is found that for 90 % of the clusters the 
deviation between the two results for the flux measure- 
ment is less than 1 %. Luminosities in the source rest 
frame energy range 0.1 — 2.4 keV have then been calcu- 



lated within XSPEC by adjusting the normalization to 
reproduce the initial flux measurements. 

The improvements of the flux determination per- 
formed here compared to the input catalogs in general 
are now summarized. 1) Due to the use of a high frac- 
tion of pointed observations the photon statistics is on 
average much better, e.g., for the 33 clusters contained 
in REFLEX and HIFLUGCS one finds a mean of 841 
and 19580 source photons, respectively. Consequently 
the cluster emission has been traced out to larger radii 
for HIFLUGCS. 2) The higher photon statistics has al- 
lowed a proper exclusion of contaminating point sources 
(stars, AGN, etc.) and substructure, and the separation 
of double clusters. 3) An iterative background correction 
has been performed. 4) A measured X-ray temperature 
has been used for the flux calculation in most cases. 

Simulations have shown that even for the HIFLUGCS 
clusters with the lowest number of photons the deter- 
mined flux shows no significant trend with redshift in 
the relevant redshift range (Ikebe et al. 2001). 
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Table 3 
Cluster properties 



Cluster 
(1) 


R.A. 

(2) 


Dec. 
(3) 




z 

(4) 


(5) 


Cx 
(6) 


A 

(7) 


»'x 
(8) 


/x 
(9) 


ix 
(10) 


(11) 


Obs 
(12) 


Ref 
(13) 


A0085 


10, 


,4632 


-9, 


.3054 


0, 


,0556 


3, 


.58 


3. 


,488 


0, 


,6 


2, 


.13 


7. 


,429 


9, 


.789 


24, 


,448 


P 


2 


A0119 


14, 


,0649 


-1, 


.2489 


0, 


,0440 


3, 


.10 


1. 


,931 


0, 


,9 


2, 


.68 


4. 


,054 


3, 


.354 


7, 


,475 


P 


1 


A0133 


15, 


,6736 


—21, 


.8806 


0, 


,0569 


1, 


.60 


1. 


,058 


0, 


,8 


1, 


.52 


2, 


,121 


2, 


.944 


5, 


,389 


P 


3 


NGC507 


20, 


,9106 


33, 


.2553 


0, 


,0165 


5, 


.25 


1. 


,093 


1, 


,3 


0, 


.88 


2, 


,112 


0, 


.247 


0, 


,326 


P 


5 


A0262 


28, 


,1953 


36, 


.1528 


0, 


,0161 


5, 


.52 


4. 


,366 


3, 


,8 


1, 


.48 


9, 


,348 


1, 


.040 


1, 


,533 


R 


6 


A0400 


44, 


,4152 


6, 


.0170 


0, 


,0240 


9, 


.38 


1. 


,146 


1, 


,1 


1, 


.85 


2, 


,778 


0, 


.686 


1, 


,033 


P 


8 


A0399 


44, 


,4684 


13, 


.0462 


0, 


,0715 


10, 


.58 


1. 


,306 


5, 


,4 


3, 


.18 


3, 


,249 


7, 


.070 


17, 


,803 


R 


6 


A0401 


44, 


,7384 


13, 


.5796 


0, 


,0748 


10, 


.19 


2. 


,104 


1, 


,1 


3, 


.81 


5, 


,281 


12, 


.553 


34, 


,073 


P 


6 


A3112 


49, 


,4912 


-44, 


.2367 


0, 


,0750 


2, 


.53 


1. 


,502 


1, 


,1 


2, 


.18 


3, 


,103 


7, 


.456 


16, 


,128 


P 


1 


FORNAX 


54, 


,6686 


—35, 


.3103 


0, 


,0046 


1, 


.45 


5. 


,324 


5, 


,6 


0, 


.53 


9, 


,020 


0, 


.082 


0, 


,107 


P+R 


4 


2A0335 


54, 


,6690 


9, 


.9713 


0, 


,0349 


18, 


.64 


3. 


,028 


0, 


,8 


1, 


.54 


9, 


,162 


4, 


.789 


7, 


,918 


P 


10 


inZw54 


55, 


,3225 


15, 


.4076 


0, 


,0311 


16 


.68 


0. 


,708 


7, 


,7 


1, 


.27 


2, 


,001 


0, 


.831 


1, 


,226 


R 


11 


A3158 


55, 


,7282 


—53, 


.6301 


0, 


,0590 


1, 


.06 


1. 


,909 


1, 


,5 


1, 


.94 


3, 


,794 


5, 


.638 


12, 


,779 


P 


1 


A0478 


63, 


,3554 


10, 


.4661 


0, 


,0900 


15, 


.27 


1. 


,827 


0, 


,6 


3, 


.12 


5, 


,151 


17, 


.690 


49, 


,335 


P 


6 


NGC1550 


64, 


,9066 


2, 


.4151 


0, 


,0123 


11, 


.59 


1. 


,979 


5, 


,4 


0, 


.71 


4, 


,632 


0, 


.302 


0, 


,407 


R 


13 


EXO0422 


66, 


,4637 


-8, 


.5581 


0, 


,0390 


6 


.40 


1. 


,390 


6, 


,2 


1, 


.32 


3, 


,085 


2, 


.015 


3, 


,283 


R 


10 


A3266 


67, 


,8410 


—61, 


.4403 


0, 


,0594 


1, 


.48 


2. 


,879 


0, 


,7 


2, 


.99 


5, 


,807 


8, 


.718 


23, 


,663 


P 


4 


A0496 


68, 


,4091 


—13, 


.2605 


0, 


,0328 


5, 


.68 


3. 


,724 


0, 


,7 


1, 


.78 


8, 


,326 


3, 


.837 


7, 


,306 


P 


8 


A3376 


90, 


,4835 


—39, 


.9741 


0, 


,0455 


5, 


.01 


1. 


,115 


1, 


,4 


2, 


.86 


2, 


,450 


2, 


.174 


4, 


,077 


P 


4 


A3391 


96, 


,5925 


—53, 


.6938 


0, 


,0531 


5, 


.42 


0. 


,999 


1, 


,9 


1, 


.98 


2, 


,225 


2, 


.681 


5, 


,857 


P 


4 


A3395S 


96, 


,6920 


-54, 


.5453 


0, 


,0498 


8, 


.49 


0. 


,836 


3, 


,8 


1, 


.45 


2, 


,009 


2, 


.131 


4, 


,471 


P 


4 


A0576 


110, 


,3571 


55, 


.7639 


0, 


,0381 


5, 


.69 


1. 


,374 


6, 


,8 


2, 


.32 


3, 


,010 


1, 


.872 


3, 


,518 


R 


6 


A0754 


137, 


,3338 


-9, 


.6797 


0, 


,0528 


4, 


.59 


1. 


,537 


1, 


,6 


1, 


.91 


3, 


,366 


3, 


.990 


11, 


,967 


P 


6 


HYDRA-A 


139, 


,5239 


-12, 


.0942 


0, 


,0538 


4, 


.86 


2. 


,179 


0, 


,6 


1, 


.66 


4, 


,776 


5, 


.930 


11, 


,520 


P 


13 


A1060 


159, 


,1784 


-27, 


.5212 


0, 


,0114 


4, 


.92 


4. 


,653 


3, 


,3 





.95 


9, 


,951 


0, 


.554 


0, 


,945 


R 


6 


A1367 


176, 


,1903 


19, 


.7030 


0, 


,0216 


2, 


.55 


2. 


,947 


0, 


,8 


1, 


.55 


6, 


,051 


1, 


.206 


2, 


,140 


P 


8 


MKW4 


181, 


,1124 


1, 


.8962 


0, 


,0200 


1, 


.86 


1. 


,173 


1, 


,7 


1, 


.23 


2, 


,268 


0, 


.390 


0, 


,543 


P 


10 


ZwC11215 


184, 


,4220 


3, 


.6604 


0, 


,0750 


1, 


.64 


1. 


,081 


1, 


,3 


2, 


.55 


2, 


,183 


5, 


.240 


11, 


,656 


P 


19 


NGC4636 


190, 


,7084 


2, 


.6880 


0, 


,0037 


1, 


.75 


3. 


,102 


7, 


,2 


0, 


.39 


4, 


,085 


0, 


.023 


0, 


,027 


R 


13 


A3526 


192, 


,1995 


-41, 


.3087 


0, 


,0103 


8, 


.25 


11. 


,655 


2, 


,2 


1, 


.64 


27. 


,189 


1, 


.241 


2, 


,238 


R 


15 


A1644 


194, 


,2900 


-17, 


.4029 


0, 


,0474 


5, 


.33 


1. 


,853 


5, 


,1 


1, 


.85 


4. 


,030 


3, 


.876 


7, 


,882 


R 


8 


A1650 


194, 


,6712 


-1, 


.7572 


0, 


,0845 


1, 


.54 


1. 


,218 


6, 


,6 


3, 


.17 


2, 


,405 


7, 


.308 


17, 


,955 


R 


6 


A1651 


194, 


,8419 


-4, 


.1947 


0, 


,0860 


1, 


.71 


1. 


,254 


1, 


,2 


2, 


.03 


2, 


,539 


8, 


.000 


18, 


,692 


P 


22 


COMA 


194, 


,9468 


27, 


.9388 


0, 


,0232 


0, 


.89 


17. 


,721 


1, 


,4 


4, 


.04 


34. 


,438 


7, 


.917 


22, 


,048 


R 


8 


NGC5044 


198, 


,8530 


-16, 


.3879 


0, 


,0090 


4, 


.91 


3. 


,163 


0, 


,5 


0, 


.56 


5, 


,514 


0, 


.193 


0, 


,246 


P 


24 


A1736 


201, 


,7238 


-27, 


.1765 


0, 


,0461 


5, 


.36 


1. 


,631 


6, 


,3 


2, 


.47 


3, 


,537 


3, 


.223 


5, 


,682 


R 


25 


A3558 


201, 


,9921 


-31, 


.5017 


0, 


,0480 


3, 


.63 


3. 


,158 


0, 


,5 


2, 


.11 


6, 


,720 


6, 


.615 


14, 


,600 


P 


1 


A3562 


203, 


,3984 


-31, 


,6678 


0, 


,0499 


3, 


.91 


1, 


,367 


0, 


,9 


2, 


.01 


2, 


,928 


3, 


.117 


6, 


,647 


P 


4 


A3571 


206, 


,8692 


-32 


,8553 


0, 


,0397 


3 


.93 


5, 


,626 


0, 


,7 


2 


.35 


12, 


,089 


8, 


.132 


20, 


,310 


P 


21 


A1795 


207, 


,2201 


26, 


,5944 


0, 


,0616 


1 


.20 


3, 


132 


0, 


,3 


2 


,14 


6, 


,270 


10, 


.124 


27, 


,106 


P 


6 


A3581 


211, 


,8852 


-27, 


,0153 


0, 


,0214 


4, 


.26 


1, 


,603 


3, 


,2 





,64 


3, 


,337 


0, 


,657 


0, 


,926 


P 


28 


MKW8 


220, 


,1596 


3, 


,4717 


0, 


,0270 


2, 


.60 


1, 


255 


8, 


A 


1 


,90 


2, 


.525 


0, 


,789 


1, 


,355 


R 


29 


A2029 


227, 


,7331 


5, 


,7450 


0, 


,0767 


3, 


.07 


3, 


294 


0, 


,6 


2 


,78 


6, 


,938 


17, 


,313 


50, 


,583 


P 


6 


A2052 


229, 


,1846 


7, 


,0211 


0, 


,0348 


2, 


.90 


2, 


,279 


1, 


,0 


1, 


,14 


4, 


,713 


2, 


,449 


4, 


,061 


P 


6 


MKW3S 


230, 


,4643 


7, 


,7059 


0, 


,0450 


3, 


.15 


1, 


,578 


1, 


,0 


1, 


,39 


3, 


,299 


2, 


,865 


5, 


,180 


P 


10 


A2065 


230, 


,6096 


27, 


,7120 


0, 


,0721 


2, 


.84 


1, 


,227 


6, 


,1 


3, 


,09 


2, 


,505 


5, 


,560 


12, 


,271 


R 


6 


A2063 


230, 


,7734 


8, 


,6112 


0, 


,0354 


2, 


.92 


2, 


,038 


1, 


,3 


2, 


,13 


4, 


,232 


2, 


,272 


4, 


,099 


P 


8 


A2142 


239, 


,5824 


27, 


,2336 


0, 


,0899 


4, 


.05 


2, 


,888 


0, 


,9 


3, 


,09 


6, 


,241 


21, 


,345 


64, 


,760 


P 


6 


A2147 


240, 


,5628 


15, 


,9586 


0, 


,0351 


3, 


.29 


2, 


,623 


3, 


,2 


1, 


.87 


5, 


,522 


2, 


,919 


6, 


,067 


P 


8 


A2163 


243, 


,9433 


-6, 


,1436 


0, 


,2010 


12, 


.27 


0, 


,773 


1, 


,5 


3, 


.15 


2, 


,039 


34, 


,128 


123, 


,200 


P 


31 


A2199 


247, 


,1586 


39, 


,5477 


0, 


,0302 


0, 


.84 


5, 


,535 


1, 


,8 


2, 


.37 


10, 


,642 


4, 


,165 


7, 


,904 


R 


8 


A2204 


248, 


,1962 


5, 


,5733 


0, 


,1523 


5, 


.94 


1, 


,211 


1, 


,6 


3, 


.29 


2, 


,750 


26, 


,938 


68, 


,989 


P 


6 


A2244 


255, 


,6749 


34, 


,0578 


0, 


,0970 


2, 


.07 


1, 


,034 


2, 


,1 


2, 


.64 


2, 


,122 


8, 


,468 


21, 


,498 


P 


6 


A2256 


255, 


,9884 


78, 


,6481 


0, 


,0601 


4, 


.02 


2, 


,811 


1, 


,4 


3, 


.09 


6, 


,054 


9, 


,322 


22, 


,713 


P 


6 


A2255 


2.58, 


,1916 


64, 


,0640 


0, 


,0800 


2, 


.51 


0, 


,976 


1, 


,2 


3, 


.22 


2, 


,022 


5, 


,506 


13, 


,718 


P 


6 


A3667 


303, 


,1362 


-56, 


,8419 


0, 


,0560 


4, 


.59 


3, 


,293 


0, 


,7 


2, 


.81 


7, 


,201 


9, 


,624 


24, 


,233 


P 


1 


SllOl 


348, 


,4941 


-42, 


,7268 


0, 


,0580 


1, 


.85 


1, 


,237 


0, 


,9 


1, 


.64 


2, 


,485 


3, 


,597 


5, 


,939 


P 


35 


A2589 


350, 


,9868 


16, 


,7753 


0, 


,0416 


4, 


.39 


1, 


,200 


1, 


,3 


1, 


.46 


2, 


,591 


1, 


,924 


3, 


,479 


P 


37 


A2597 


351, 


,3318 


-12, 


,1246 


0, 


,0852 


2, 


.50 


1, 


,074 


1, 


,2 


1, 


.43 


2, 


,213 


6, 


,882 


13, 


,526 


P 


6 


A2634 


354, 


,6201 


27, 


,0269 


0, 


,0312 


5, 


.17 


1, 


,096 


1, 


,6 


1, 


.79 


2, 


,415 


1, 


,008 


1, 


,822 


P 


6 
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3.2. Mass Determination 

A parametric description of the cluster gas density 
profile has been derived using the standard (3 model (e.g., 
Cavaliere & Fusco-Femiano 1976; Gorenstein et al. 1978; 
Jones & Forman 1984). Assuming spherical symmetry 
the model 



jD2 



-3/3+1 



B 



(1) 



is fitted to the measured surface brightness profile (ring 
count rates per ring area), where R denotes the projected 
distance from the cluster center. This yields values for 
the core radius, Tc, the /3 parameter, and the normaliza- 
tion, 5*0 (and also a fitted value for the background sur- 
face brightness, B, since the fit is performed on the non 
background subtracted data). The fit values have been 
used to construct the radial gas density distribution 



r) = /oo 1 



(2) 



The gravitational cluster mass, Mtot, has been deter- 
mined assuming the intracluster gas to be in hydrostatic 

equilibrium and isothermal. Using (2) and the ideal gas 
equation under these assumptions leads to 



Mtot(< r) = 



/impG 



1 



(3) 



where jj, (= 0.61) represents the mean molecular weight, 
nip the proton mass, and G the gravitational constant. 
Combined iV-body/hydrodynamic cluster simulations 
have shown that this method generally gives unbiased 
results with an uncertainty of 14-29% (e.g., Schindler 
1996a; Evrard et al. 1996). Currently there are contra- 
dictory measurements concerning the general presence 
of gas temperature gradients in clusters (e.g., Fukazawa 
1997; Markevitch et al. 1998; Irwin et al. 1999; White 
2000; Irwin & Bregman 2000). If there is a systematic 
trend in the sense that the gas temperature decreases 
with increasing radius in the outer cluster parts similar 
to that found by Markevitch et al. then the isothermal 
assiimption leads to an overestimation of the cluster mass 
of about 30% at about 6 core radii (Markevitch et al. 
1998). Finoguenov et al. (2001) determined masses by 
employing the assTimption of isothermality and also using 
measured cluster gas temperature profiles. A comparison 
for 38 clusters included in their sample indicates that the 
latter masses are on average a factor of 0.80 smaller than 
the isothermal masses within r2oo (this radius is defined 
in the next paragraph) . Until the final answer on this is- 
sue is given by XMM- Newton"^, we retain the isothermal 



^First results indicate that apart from the central regions the in- 
tracluster gas is isothermal out to at least l/2r2oo (e.g., Reiprich 
2001b; M. Arnaud 2001, private communication). 



assumption. The influence of a possible overestimation 
of the cluster mass on the determination of cosmological 
parameters is investigated in Sect. 5.3.2. 

Having determined the integrated mass as a function 
of radius, a physically meaningful fiducial radius for the 
mass measurement has to be defined. The radii com- 
monly used are either the Abell radius, r200; or rsoo- 
The Abell radius is fixed at ta = 3 h^Q Mpc. The radius 
''200 ('"soo) is the radius within which the mean gravi- 
tational mass density (ptot) = 200 (500) pc- The critical 
cosmic matter density is defined as pc = 3H'^/{8ttG), 
where ff^ ^ E{zf and E{z) = [17^(1 + zf + 
rik(l + + r^A]^/^. It has been shown that a correc- 
tion for redshift is not necessary for the nearby clusters 
included in HIFLUGCS (Finoguenov et al. 2001) and 
we use the zero redshift value for all calculations, i.e. 
Pc = 4.6975 X 10~^°gcm~^, unless noted otherwise. Nev- 
ertheless the influence of this approximation is tested in 
Sect. 5.3.2 for the model (fim = l,f^A = 0), where evo- 
lution is strong. 

In order to treat clusters of different size in a homo- 
geneous way we determine the cluster mass at a char- 
acteristic density but also give the mass determined for- 
mally at a flxed radius for comparison. Spherical collapse 
models predict a cluster virial density (pvir) ~ 178 pc for 
(rim = 1, = 0), so a pragmatic approximation to the 
virial mass is to use rgoo as the outer boundary. Sim- 
ulations performed by Evrard et al. (1996) have shown, 
however, that isothermal X-ray mass measurements may 
be biased towards high masses for r > rsoo. Further- 
more for most of the clusters in HIFLUGCS (86%) up 
to rsoo no extrapolation outside the significantly detected 
clustcir emission is necessary, i.e. rsoo < rx, whereas 
the fraction is lower for r2oo (25%) and rA (17%). In 
summary the most accurate results are expected for 
Mtot{< fzao) = Afsoo, but for a comparison to pre- 
dicted mass functions M200 is the more appropriate value 
(Sect. 5.3.2). Results for all determined masses and their 
corresponding radii are given in Tab. 4. Masses for the 
cluster gas will be given in a subsequent paper. 

A major source of uncertainty comes from the tem- 
perature measurements. However, this (statistical) error 
is less than 5 % for one third of the clusters, therefore 
also other sources of error have to be taken into account, 
in particular one cannot neglect the uncertainties of the 
fit parameter values when assessing the statistical errors 
of the mass measurements. Therefore mass errors have 
been calculated by varying the fit parameter values, /? 
and Vc , along their 68 % confidence level error ellipse and 
using the upper and lower bound of the quoted temper- 
ature ranges. The statistical mass error range has then 
been deflned between the maximum and minimum mass. 
Note that a simple error propagation applied to (3) would 
underestimate the uncertainty of M200 and M500, since 
r2oo and rsoo also depend on Tgas, P, and (weakly) rc- 
The individual mass errors have been used in subsequent 
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Table 3 — Continued 



Cluster 


R.A. 


Dec. 


z 


riH 


Cx 


A 




/x 


ix 


LboI 


Obs 


Ref 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


A2657 


356.2334 


9.1952 


0.0404 


5.27 


1.148 


0.9 


1.52 


2.535 


1.771 


3.202 


P 


8 


A4038 


356.9322 


-28.1415 


0.0283 


1.55 


2.854 


1.3 


1.35 


5.694 


1.956 


3.295 


P 


4 


A4059 


359.2541 


-34.7591 


0.0460 


1.10 


1.599 


1.3 


1.72 


3.170 


2.872 


5.645 


P 


36 






Clusters from the extended 


sample not included in 


HI FLU CCS. 








A2734 


2.8389 


-28.8539 


0.0620 


1.84 


0.710 


2.5 


1.74 


1.434 


2.365 


4.357 


P 


1 


A2877 


17.4796 


-45.9225 


0.0241 


2.10 


0.801 


1.2 


1.06 
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0.714 


P 


4 
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0.0147 


5.25 


0.313 


2.5 


0.30 


0.479 


0.045 


0.051 


P 


5 
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0.0172 


9.21 


7.007 


2.0 


1.58 
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3.882 


R 


7 
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40.723 


0.8 


3.30 
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16.286 
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R 


9 


S405 
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7.65 


0.781 


8.2 
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5.574 


R 


12 
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5.6 


1.61 
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2.242 


4.996 


R 


10 
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79.1560 


6.4421 
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1.221 


1.3 


1.37 
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P 
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0.788 
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R 


4 
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5.4 


0.73 


0.234 
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0.240 


P 


15 
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-25.4692 


0.0410 


1.88 


0.771 


1.8 


2.12 


1.551 
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1.870 


P 


15 
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0.0498 


5.42 
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3.9 


1.37 
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3.461 


P 


4 
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0.0340 
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0.94 
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R 


16 
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2.44 
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70.604 


P 


17 
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5.14 
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P 
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157.5151 


-35.3093 


0.0116 
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R 


18 
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P 


6 
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1.59 
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1.0 


0.27 
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P 


15 
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P 
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6.00 
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1.55 
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2.317 


P 


21 
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0.0539 


5.96 


0.797 


1.8 


1.64 
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4.483 


P 


21 
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-1.3408 


0.1840 


1.80 


0.712 


1.1 


2.36 


1.454 


20.605 


60.707 


P 


23 
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0.0495 


3.92 


0.714 


2.5 


2.00 
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1.601 


2.701 


P 


26 
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0.0757 


1.00 


0.654 


1.8 


2.02 
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P 


27 
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207.3408 


28.1038 


0.0748 


1.18 


0.610 


7.9 


1.98 


1.183 


2.840 


5.337 


R 


28 
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0.1712 


0.97 


0.729 


1.4 


2.35 


1.454 


17.813 


56.533 


P 


6 


NGC5813 
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1.6981 


0.0064 


4.19 


0.976 


6.7 


0.17 


1.447 


0.025 


0.029 


R 


13 


NGC5846 


226.6253 


1.6089 


0.0061 


4.25 


0.569 


2.3 


0.21 


0.851 


0.014 


0.016 


P 


13 
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0.0369 


3.36 
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1.9 


1.46 
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0.917 
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P 


8 
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20.83 
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2.20 
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R 


30 
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0.0510 


12.29 
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0.7 


2.54 
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P 


32 
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258.1115 


-23.3634 


0.0280 


20.14 


11.642 


2.0 


2.29 


35.749 


11.953 


37.391 


R 


33 
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266.0623 


32.9893 


0.0757 


3.56 


0.889 


4.4 


1.83 


1.850 


4.529 


9.727 


R 


34 
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43.9484 


0.0564 


8.77 


5.029 


1.0 


3.57 


12.202 


16.508 


47.286 


P 


6 


A3695 


308.6991 


-35.8135 


0.0890 


3.56 


0.836 


9.2 


2.58 


1.739 


5.882 
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R 


1 
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2.5564 


0.0494 


6.63 


0.841 


7.3 


2.20 


1.884 


1.969 


3.445 


R 


5 


A3822 


328.5438 


-57.8668 


0.0760 


2.12 


0.964 


7.3 


3.18 


1.926 


4.758 


9.877 


R 


1 


A3827 


330.4869 


-59.9641 


0.0980 


2.84 


0.953 


5.8 


1.78 


1.955 


7.963 


20.188 


R 


1 


A3888 


338.6255 


-37.7343 


0.1510 


1.20 


0.546 


2.4 


1.52 


1.096 


10.512 


30.183 


P 


23 


A3921 


342.5019 


-64.4286 


0.0936 


2.80 


0.626 


1.7 


2.43 


1.308 


4.882 


11.023 


P 


12 


HCG94 


349.3041 


18.7060 


0.0417 


4.55 


0.820 


1.0 


2.09 


1.775 


1.324 


2.319 


P 


36 


RXJ2344 


356.0723 


-4.3776 


0.0786 


3.54 


0.653 


1.4 


1.61 


1.385 


3.661 


7.465 


P 
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Note. — Column (1) lists the cluster name. Names have been truncated to at most eight characters to preserve the 
compactness of the tabic. Columns (2) and (3) give the equatorial coordinates of the cluster center used for the regional 
selection for the epoch ,12000 in decimal degrees. Column (4) gives the heliocentric cluster redshift. Column (5) lists the 
column density of neutral galactic hydrogen in units of lO'^" atoms cm^^. Column (6) gives the count rate in the channel 
range 52—201 which corresponds to about (the energy resolution of the PSPC is limited) the energy range 0.5 — 2.0 keV in 
units of ctss~^. Column (7) lists the relative 1-a Poissonian error of the count rate, the flux, and the luminosity in percent. 
Column (8) gives the significance radius in h^^ Mpc. Column (9) lists the flux in the energy range 0.1 — 2.4 keV in units of 
10~^^ ergss"^ cm~^. Column (10) gives the luminosity in the energy range 0.1 — 2.4 keV in units of /i^q^ 10^* ergs s~^. Column 
(11) gives the bolometric luminosity (energy range 0.01 — 40keV) in units of h^Q 10*'* ergss~*. Column (12) indicates whether 
a RASS (R) or a pointed (P) ROSAT PSPC observation has been used. Column (13) lists the code for the redshift reference 
decoded below. 

References. — (1) Katgcrt ct al. 1996. (2) Mazure et al. 1996. (3) Median of 9 galaxy redshifts compiled from Lauberts & 
Valentijn 1989; Merrifield & Kent 1991; Lovcday et al. 1996; Way ct al. 1998. (4) AbcU ct al. 1989. (5) Huchra et al. 1999. (6) 
Struble & Rood 1987. (7) dell'Antonio ct al. 1994. (8) Zabludoff ct al. 1993. (9) Poulain ct al. 1992. (10) NED Team 1992. 
(11) Bohringer et al. 2000. (12) Dc Grandi ct al. 1999. (13) de Vaucoulcurs ct al. 1991. (14) Zabludoff et al. 1990. (15) den 
Hartog & Katgert 1996. (16) Michel & Huchra 1988. (17) Yan & Cohen 1995. (18) Garcia 1995. (19) EbeUng ct al. 1998. 
(20) Median of 8 galaxy redshifts compiled from de Vaucoulcurs et al. 1991; Quintana et al. 1995; Katgcrt et al. 1998. (21) 
Vcttolani et al. 1990. (22) Allen ct al. 1992. (23) Teague et al. 1990. (24) da Costa et al. 1998. (25) Dressier & Shcctman 
1988. (26) Melnick & Moles 1987. (27) Median of 13 galaxy redshifts compiled from Kirshner et al. 1983; Zabludoff et al. 1990; 
NED Team 1992; Davoust & Considere 1995; Oegerlc et al. 1995. (28) Postman et al. 1992. (29) Andersen & Owen 1994. (30) 



calculations, unless noted otherwise.^ A mean statistical 
error of 23% for clusters included in HIFLUGCS and 
a mean error of 27% for the extended sample has been 
found. 



®In log space errors are transformed as A log a; = log(e) — 
x~)/{2x), where a;+ and x~ denote the upper and lower boundary 
of the quantity's error range, respectively. 



11 



Table 4 
Cluster properties 
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4. Results 

In this Section it is shown that a tight correlation ex- 
ists between the gravitational cluster mass and the X-ray 
luminosity. This ensures that HIFLUGCS is essentially 
selected by cluster mass. In the second part of this Sec- 
tion the cluster mass function is prcsc^ntcd, including the 
proper treatment of the scatter in the Lx--^tot relation. 

4.1. Mass— Luminosity Relation 

Since the aim is the construction of a mass function 
from a flux-limited sample it is now important to test for 
a correlation between X-ray luminosity and gravitational 
mass. In Fig. 6 Lx, given in the ROSAT energy band, is 
plotted as a function of M200, showing clearly the exis- 
tence of a tight (linear Pearson correlation coefficient = 
0.92) correlation, as expected. 
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Fig. 6. — Gravitational mass-X-ray luminosity relation 
(solid line) for the extended sample of 106 galaxy clus- 
ters. The dashed line gives the best fit relation for the 
63 clusters included in HIFLUGCS (filled circles only). 
The bisector fit results are shown. One-cr statistical error 
bars are plotted for both axes, however, only the mass 
errors are larger than the symbol sizes. 

To quantify the mass-luminosity relation, a linear re- 
gression fit in log-log space has been performed. The 

method used allows for intrinsic scatter and errors in 
both variables (Akritas & Bershady 1996). Tables 7-11 
in the appendix give the results for different fit methods, 
where minimization has been performed in vertical, hor- 
izontal, and orthogonal direction, and the bisector result 
is given, which bisects the best fit res\.ilts of vertical and 
horizontal minimization. The fits have been performed 
using the form 



log 



Lx(0.1 -2.4keV)\ 
/i^o^ 1040 ergs s-i / 



= ^ + a log 



200 



(4) 



We find, as noted in general by previous authors (e.g., 
Isobe et al. 1990), that the chosen fitting method has 



a significant influence on the best fit parameter values.^ 
In this work the appropriate relation for the application 
under consideration is always indicated. 

The difference between the fit results for 63 and 106 
clusters may indicate a scale dependence of the Lx- 
Mtot relation, since the difference is slightly larger than 
the uncertainty evaluated with the bootstrap method. 
The small mimbcr of low luminosity clusters in HI- 
FLUGCS compared to the extended sample may be re- 
sponsible for the less steep relation obtained using the 
HIFLUGCS clusters only. Note that only two out of 
the six clusters with Lx < h^Q 10^^ ergss"^ are included 
in HIFLUGCS. To reliably detect any deviations from 
the power law shape of the Lx-Mtot relation, however, 
more clusters with Mtot < 10^^ h^Q Mq (and possibly 
Mtot > 3 X 10^^ h'^Q Mq) need to be sampled. Such work 
is in progress. As will be seen later, in the procedure 
used here for the comparison of observed and predicted 
mass functions the precise shape of the Lx-Mtot relation 
is not important. 

When constructing the mass function the overall 
(measurement plus intrinsic) scatter in the Lx^Aftot re- 
lation may become important (Sect. 4.2). After verifying 
that the scatter is approximately Gaussian in log space 
the scatter has been measured as given in Tab. 12 in 
the Appendix. The scatter in log(Lx(0.1 - 2.4 keV)), 
log(M2oo), and orthogonal to the best fit line is given by 
o-iogLx) criogMtot) and aiogL/M, respectively. 

4.2. Mass Function 

The commonly used definition of the galaxy cluster 
mass function is analogous to the definition of the lumi- 
nosity function (e.g., Schechter 1976): the mass func- 
tion, (/)(M), denotes the number of clusters, N, per 
unit comoving volume, dV, per unit mass in the in- 
terval [M,M + dM], i.e. (j){M) = N{M)/{dVdM) = 
dn{M)/dM. Assuming constant density the classical 
Vmax estimator (e.g., Schmidt 1968; Felten 1976; Binggeli 
et al. 1988) can be used for estimation of luminosity func- 
tions, i.e. (j>{L) ~ 1/AL J2iLi ■ ■ ^max is the maxi- 
mum comoving volume within which a cluster with given 
luminosity for a given survey flux limit and sky cover- 
age could have been detected. As mentioned in Sect. 2 
the HIFLUGCS flux limit is constant over 99% of the 
covered area, which simplifies the calculation of V^ax- 

In the previous Section it has been shown that the 
X-ray luminosity is closely correlated with cluster mass. 
Therefore the V^ax estimator can also be applied to es- 
timate the mass function; Vmax then being a function of 
mass. Wc employ different methods to correct for the 



^This also implies that for a proper comparison of relations which 
have been quantified by many different authors, e.g., the ix" Tx 
relation, one and the same fitting statistic ought to be used (e.g., 
Wu et al. 1999). 
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Note. — Column (1) lists the cluster name. Golumn (2) gives rhe /3 parameter value and the corresponding 68% c.l. statistical 
uncertainty for t-wo interesting parameters. Column (3) gives the core radius in h'^Q kpc and the corresponding uncertainty. 
Column (4) lists the X-ray temperature along with its error. For some references the temperature uncertainty is quoted at the 
90% confidence level and therefore represents a conser-vative error estimate. Columns (5) and (7) give Msoo and M200 and 



scatter present in the Lx^^tot relation. If VmaxiLx) 
is used instead of Kiax(^^tot) = Kiax(-^(^^tot)), where 
L(Mtot) is the luminosity estimated from the I/x--^tot 
relation using the determined cluster mass Mtot, the scat- 
ter is automatically taken into account. This method 
has been widely used in the construction of X-ray tem- 
perature functions, recently, e.g., by Henry (2000). If 
Knax(Mtot) is used and the utilized Lx-Mtot relation is 
assumed to be the 'true' relation then the scatter in this 
relation has to be taken into account explicitly. There- 
fore following the method employed for the temperature 
function by Markevitch (1998) and Ikebe et al. (2001) 
the mass function may also be estimated by determin- 
ing V^g^^{Mtot), where the measured scatter in logLx is 
included. Specifically we use 

/oo 
-OO 

X exp dlogL , (5) 

where A and a are the best fit parameter values taken 
from the appropriate Lx-Mtot relation of the form (4) 
and <T\ogLx is the corresponding measured standard de- 
viation in logZ/x given in Tab. 12. However, we can 
also use the measured Lx-Mtot relation directly, i.e. 
Vmax(-^tot), taking advantage of the fact that in our flux- 
limited sample there are fewer low luminosity clusters for 
a given mass than high luminosity ones, which results in 
a slightly increased normalization of the relation. There- 
fore using this relation directly, the effect of the scatter 
and the resulting bias towards higher luminosity clusters 
is already included and thus directly accounted for. 

The drawback of using Vmax{Lx) is that a small num- 
ber of clusters per mass bin possibly does not represent 
the true scatter well. To minimize this effect we use at 
least ten clusters per mass bin. The drawback of using 
V^^^{Mtot) or K„ax(Mtot), as noted, e.g., by Markevitch 
(1998) and Eke et al. (1998), is the reliance on the valid- 
ity of the measured relation over the entire mass range. 
The first method and the method that accounts for the 
scatter explicitly (eq. 5) have been tested by using Monte 
Carlo simulations for a precisely known Lx-Tx relation 
and scatter and have been shown to give accurate esti- 
mates of <j){T) for a large number of clusters in the study 
of the HIFLUGCS temperature function by Ikebe et al. 
(2001). 

In Fig. 7 HIFLUGCS mass functions arc shown. As 
expected the method employing Vmax(ix) prompts a 
mass function exhibiting a larger scatter, because in this 
case the scatter is accounted for by the actual scatter of 
the ten or eleven clusters in each mass bin. For com- 
parison the two extreme mass functions calculated using 
Vmax(Aftot) are shown. Extreme is meant in the sense 
of using the steepest (A) and shallowest (B) Lx -Mtot 
relation for the HIFLUGCS sample, i.e., (M | L) with 



a = 1.538 and (L | M) with a = 1.310 (Tab. 8). At the 

low mass end (A) predicts a lower luminosity for a given 
mass than (B) resulting in a smaller V^nax and therefore 
a higher dn/dM. At the high mass side the effect is 
opposite resulting in a lower dn/dM for (A). The differ- 
ences of these mass functions to the mass function calcu- 
lated using Vma,x{Lx) can be understood in a similar way 
and are caused partly by the indication of a deviation 
from a power law shape of the Lx-Mtot relation. Using 
V,*ax(-^tot) results in similar mass functions as shown for 
the open symbols in Fig. 7 but the points lie systemat- 
ically lower because the scatter is accounted for twice. 
For the comparison of the observational mass function to 
mass functions predicted by certain cosmological models 
Knax(i'x) is used because it is independent of the pre- 
cise shape of the Lx Mtot relation and also because Lx 
has a much smaller measurement uncertainty than Mtot ■ 
The influence of the choice of the Vmax calculation on the 
estimation of cosmological parameters is investigated in 
Sect. 5.3.2. 
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Fig. 7. — Gravitational mass functions for HIFLUGCS. 

The mass function plotted with filled circles has been 
determined using T4nax(-^x), the ones with open symbols 
using Viiiax(A'/tot), where triangles correspond to the re- 
lation (L I M) and squares to (M | L) both for the flux- 
limited sample (see text Sect. 4.2). Vertical error bars 
correspond to the formal 1-a Poisson errors, the horizon- 
tal bars indicate the mass intervals covered. Each bin 
contains 10 clusters, apart from the highest mass bin, 
which contains 11 clusters. The highest and lowest mass 
clusters have been used to calculate the highest and low- 
est mass intervals. 



5. Discussion 

A precise determination of distribution functions re- 
quires a high sample completeness. In Sect. 5.1 several 
completeness tests for HIFLUGCS are discussed, indi- 
cating a high completeness. The observed Lx-Mtot rela- 
tion is compared to expectations in Sect. 5.2 and possible 
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applications are indicated. The cluster mass function is 
compared to previous determinations and to predictions 
of cosmological models in Sect. 5.3. The total gravita- 
tional mass contained in galaxy clusters is compared to 
the total mass in the universe in Sect. 5.4. 

5.1. Sample Completeness 

The sample completeness is important for the accu- 
racy of the mass function. The selection criteria detailed 
in Sect. 2 are met by 63 clusters with mean redshift 
{z) = 0.05 and with two clusters having z > 0.1. The 
sample is constructed from surveys with much deeper 
flux limits and high completenesses. A possible remain- 
ing incompleteness in these surveys is likely to be present 
at low fluxes close to their flux limits, which therefore 
would not eSect HIFLUGCS. Nevertheless four com- 
pleteness tests have been performed and are described in 
this Section; they all indicate a high completeness of HI- 
FLUGCS. The logA^-logS* and Lx-z diagram are com- 
pared to expectations, the luminosity function is com- 
pared to luminosity functions of deeper surveys, and the 
V/Vmax test is performed. 
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Fig. 8. -- \ogN{> /x)-log/x diagram. Fluxes are mea- 
sured in the ROSAT energy band (0.1 - 2.4keV). The 
dashed line has a slope —1.5, expected for a uniform 
distribution of clusters in a static Euclidean universe 
('three-halves-law'), the line is normalized to produce the 
same cluster number at /x = 8 x 10~^^ ergss"^ cm~^. 

Figure 8 shows the integral number counts as a func- 
tion of X-ray flux ('log N log S") . The slope in the log A^- 
logS* diagram is very close to the value —1.5 expected 
in a static Euclidean universe for uniformly distributed 
clusters. Due to the small number of clusters (4) the de- 
viation is not significant for /x ^ 1 x 10~^° ergs s~^ cm~^. 
Since the average redshift is smallest for the highest 
fluxes large scale structure is not completely washed out 



at the high flux end, therefore the slight bump visible 
around /x 6 x 10^^^ ergs s~^ cm^^ in Fig. 8 suggests 
a deviation caused by cosmic variance. The effect of an 
expanding and finite universe on the log A'-log S - fiat- 
tening of the slope towards low fluxes is small for the 
redshift range covered by the sample. The slope consis- 
tent with —1.5 towards the flux limit therefore indicates 
a high completeness of HIFLUGCS. 

In Fig. 9 the X-ray luminosity is plotted as a function 
of redshift. The flux limit is shown as a solid line.^ One 

" 10^ " • ' • • • • • • ' • • • 
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Fig. 9. X-ray luminosity as a function of redshift. The 
flux limit is shown as a solid line. 

notes the increase in rare luminous systems with increas- 
ing redshift (volume). Because of the seeming underden- 
sity of clusters in the redshift range 0.10 < z < 0.15 a 
comparison with the expected number of clusters as de- 
rived from A-body simulations has been performed. An 
OCDM simulation, carried out for analysis of the power 
spectral densities of REFLEX clusters (Schuecker et al. 
2001b), adjusted to the HIFLUGCS survey volume in 
the southern hemisphere (roughly half of the total vol- 
ume sampled by HIFLUGCS) has been used. The sim- 
ulation yields 39 clusters while 33 HIFLUGCS clusters 
have been detected in this region. It is found that in fact 
not even one cluster with z > 0.1 is expected for this 
volume based on this simulation and the HIFLUGCS 
subsample also does not contain any cluster with a red- 
shift larger than 0.1. This is a further piece of evidence 
for the high completeness of the sample. 

In Fig. 10 the HIFLUGCS X-ray luminosity func- 
tion is compared to luminosity functions of larger surveys 
in the southern (REFLEX, Bohringer et al. 2001a) and 

^The correction i<"(Tgas, z) for converting observer rest frame lumi- 
nosities to source rest frame luminosities depends on redshift and 
source spectrum (Tgas). For source rest frame luminosities it is 
therefore not possible to plot the flux limit as one line in 2 dimen- 
sions {Lx,z), but rather as an area in 3 dimensions (Lx,z,Tgas). 
For consistency we therefore give in this 2d plot the observer rest 
frame luminosity (the correction is less than 6 % for 90 % of the 
clusters anyway). 
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northern (BCS, Ebeling et al. 1997) hemisphere. These 
surveys have much deeper flux Hmits (Sect. 2) and con- 
tain many more clusters. Very good agreement is found, 
which shows the high completeness and homogeneous se- 
lection of HIFLUGCS. 
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Fig. 10. — X-ray luminosity function for HIFLUGCS 
compared to surveys with deeper flux limits in the north- 
ern (N) and southern (S) hemisphere. Vertical error bars 
correspond to the formal 1-a Poisson errors (no cosmic 
variance), the horizontal bars indicate the luminosity in- 
tervals covered. 

The V/Vmax test (e.g.. Rowan- Robinson 1968; Schmidt 

1968; Avni & Bahcall 1980; Peacock 1999, Sect. 14.5) 
can be used to asses a possible sample incomplete- 
ness. Assuming a uniform distribution of clusters a 
value 1/2 is expected on average. For HIFLUGCS 
(V/Vmax) = 0.47 ± 0.04, which is consistent with the 
expectation and we interpret this rc!sult as a clear sign 
that HIFLUGCS covers a large enough volume for most 
of the Lx range to be representative of the local universe 
with a high sample completeness. The local nature of 
HIFLUGCS becomes obvious by noting that the result 
of the comoving V/Vma.^ test is almost identical to the 
result of the equivalent test assuming a Euclidean and 
non expanding space, i.e., ((/x//x,iim)"^/^) = 0.46. 

5.2. Mass— Luminosity Relation 

The close correlation between the X-ray luminosity 
and the gravitational mass found in Sect. 4.1 is not 
surprising. Simple self similar scaling relations predict 
Tgas oc Mtoti?^^ and Mtot oc R^^^, where i?ch is a char- 
acteristic radius, e.g., the virial radius. Combined with 
bremsstrahlung emission, LboI oc p^^T^^'^R^^ (Lx(0.1 ~ 

2.4 keV) « P^as-R?h), the relation LboI « /^..Af^'/t' 
(Lx(0.1 - 2.4 keV) oc f^^Mtot) is predicted (e.g., Per- 
renod 1980), where the gas fraction /gas = MgasMj"^ 
and LboI is the bolometric luminosity. 

Observationally from the tight correlations between X- 



ray luminosity and temperature (e.g., Markevitch 1998), 
and temperature and mass (e.g., Finoguenov et al. 2001) 
a correlation between luminosity and mass clearly is ex- 
pected. Also correlations found between X-ray luminos- 
ity and galaxy velocity dispersion (e.g.. Edge & Stewart 
1991b) and X-ray luminosity and mean shear strength 
from weak lensing studies (e.g.. Small et al. 1997) indi- 
cate a correlation between Lx and Mtot- 

The X-ray luminosity has been compared directly to 
gravitational mass estimates by Reiprich (1998), Reiprich 
& Bohringer (1999a, 2000), Schindler (1999), Jones & 
Forman (1999), Miller et al. (1999), Ettori & Fabian 
(2000), and Borgani & Guzzo (2001), where good cor- 
relations have been found in all of these studies. 
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Fig. 11. — Gravitational mass-bolometric X-ray lumi- 
nosity relation compared to predicted relations. Shown 
are: best fit relation for the extended sample (solid line), 
best fit relation determined using HIFLUGCS (triple- 
dot-dashed line), self-similar relation normalized by sim- 
ulations of Navarro et al. (1995) (dot-dashed line), pre- 
heated relation given by Evrard & Henry (1991), using a 
normalization taken from the simulations of pre-heated 
clusters by Navarro et al. (1995) (dashed line). 

In order to compare the empirical Lx A^toi relation 
with predictions a quasi bolometric luminosity, LboI, has 
been calculated in the source rest frame energy range 
0.01—40 keV (for the relevant range of cluster gas temper- 
atures at least 99 % of the flux is contained in this energy 
range). In Fig. 11 this Z/bo1-^200 relation is compared 
to predicted relations. The solid line shows the best flt 
relation for the 106 clusters in the extended sample and 
the triple-dot-dashed line shows the best fit relation de- 
termined using HIFLUGCS. Here the bisector fit results 
have been used in order to treat variables symmetrically, 
which is the appropriate method for a comparison to the- 
ory (e.g., Isobe et al. 1990). The dot-dashed line shows 
the self-similar relation (LboI oc Af^/^) normalized by the 
simulations of Navarro et al. (1995) and the dashed line 
shows the 'pre-heated' relation given by Evrard & Henry 
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(1991) (-LboI oc M^^/^), using a normalization taken from 
the simulations of pre-heated clusters by Navarro ct al. 
(1995). The idea of pre-heating is that the intracluster 
gas is not cold initially, as in the self similar case, but 
is heated by some form of non gravitational heat input, 
e.g., from supernovae or AGN, before or during cluster 
formation. Assuming the central regions of all clusters 
to have the same entropy yields the latter relationship. 
Fig. 11 shows that measured and predicted relations are 
in rough agreement, the difference between the predicted 
relations being larger than the difference to the observed 
relations. Note, however, that the X-ray luminosity is 
one of the most uncertain quantities to be derived from 
simulations. Frenk et al. (1999) recently showed in a com- 
parison of twelve different cosmological hydrodynamics 
codes that a factor of two uncertainty is a realistic es- 
timate of the current accuracy. Including gas cooling 
in simulations worsens the situation (e.g., Balogh et al. 
2001). The slopes of the observed relations are closer to 
the pre-heated relation. Observationally the effect of pre- 
heating can also result in a decrease of the gas mass frac- 
tion for low temperature systems. This has actually been 
observed for the clusters in our sample (Reiprich 1998; 
Rciprich & Bohringcr 1999b). The possibility that winds 
from, e.g., supernovae - originally invoked to explain the 
apparent low gas content of elliptical galaxies (e.g., Math- 
ews & Baker 1971; Larson 1974) prc-heat and dilute 
the central gas and thereby break the self similarity has 
been pointed out by various authors (e.g.. Kaiser 1986). 
Such a process would work most efficiently on the least 
massive clusters (e.g.. White III 1991; Mctzler & Evrard 
1997; Ponman et al. 1999). The Lx M^^t relation and 
other relevant relations between physical cluster param- 
eters of the extended sample have been discussed more 
thoroughly in this context in Reiprich (2001a). 

One may wonder about the origin of the scatter in 
the Mtot-Lx relation. Due to the use of pointed ob- 
servations for most of the clusters and the local nature 
of HIFLUGCS the statistical errors of Lx are negligi- 
ble. The logarithmic mean mass measurement uncer- 
tainty has been measured as 0.12. The overall scatter in 
log mass of the data points compared to the best fit re- 
lation is larger and has been measured as 0.21 (Tab. 12). 
This indicates a possible contribution of intrinsic scatter 
to the overall scatter in the Mtot-Lx relation. An obvious 
candidate to cause intrinsic scatter is the central excess 
emission (central surface brightness exceeding single /? 
model surface brightness) present in a number of clusters. 
This excess emission may have its physical origin either 
in cooling flows (e.g., Fabian 1994) or in the presence of 
cD galaxies (e.g., Makishima et al. 2001). A cooling flow 
analysis of the HIFLUGCS clusters is in progress and 
first results indicate that indeed clusters with a high in- 
ferred mass deposition rate lie on the high Lx side of the 
Mtot-Lx relation (Y. Chen et al., in preparation). 

In Fig. 12 the measured number of cluster member 



galaxies as taken from Abell et al. (1989) is compared 
to Lx as gravitational mass tracer. It is clearly seen 
that a selection by X-ray luminosity is much more effi- 
cient than a selection by Abell richness in terms of mass. 
Even though only the X-ray surface brightness profile 
and neither its normalization nor the X-ray luminosity 
are directly used in the X-ray mass determination via the 
hydrostatic equation, it is nonetheless reassuring that a 
similar result is obtained when Lx and richness are com- 
pared to masses estimated from optical velocity disper- 
sions (Borgani & Guzzo 2001). 
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Fig. 12. — Measured number of cluster member galax- 
ies, A^gx, as taken from Abell et al. (1989) and X-ray 
luminosity of the same (66) clusters as a function of the 
gravitational mass. Double clusters, whose components 
have been treated separately here, e.g., A3528n/s, are re- 
moved. Above the dashed line all clusters have an Abell 
richness R> 1. 

A wide range of possible applications becomes avail- 
able with the quantification of the Lx^Mtot relation and 
its scatter. For large X-ray flux-limited cluster surveys, 
where individual mass determinations are currently not 
feasible, luminosities can be directly converted to masses. 
No combination of observations, simulations and theory 
is then needed, like the frequently used approach of re- 
lating X-ray luminosities to X-ray temperatures by an 
observed relation, and converting X-ray temperatures to 
masses using a relation where the slope is taken from the- 
oretical arguments and the normalization from hydrody- 
namical simulations (e.g., Moscardini et al. 2000). The 
observational Lx-Mtot relation has first been applied di- 
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rectly in this sense in the power spectral analysis of RE- 
FLEX clusters (Schuecker et al. 2001b). An example 
of another direct application is given in Sect. 5.3.3. The 
Lx-Mtot relation may also be applied to convert theoret- 
ical or simulated mass functions to luminosity functions 
for comparison with observations of X-ray flux-limited 
samples, which is currently being performed in the inter- 
pretation of the REFLEX luminosity function. 

At this point it is important to note that even for the 
highest redshift cluster in our sample (z = 0.2) the de- 
pendence of the observational determination of Mtot and 
Lx on the chosen cosmological model is very weak. For 
instance a.t z = 0.2 the increase in the Imninosity dis- 
tance, Dl, and the diameter distance. Da, is less than 
5% going from (O^ = 1.0, I^a = 0) to {fl^ = 0.1, 17a = 
0). From (3) one finds that -Mtot(< r) oc r and therefore 
Mtot{< r) oc Da, implying an increase of Mtot by less 
than 5 % for the two models above. For Lx one has an 
increase of less than 10%. This means that the Lx -Mtot 
relation given here can be used unchanged for various cos- 
mological applications (unless redshift ranges are probed 
where evolution becomes important, in this case a model 
dependent redshift correction has to be introduced). A 
similar calculation for Knax shows that for the extreme 
case Zran-K = 0.2 the increase in Vmax is less than 14%, 
which is less than the size of the Poissonian error bars in 
Fig. 14. 

More detailed investigations on the shape of the re- 
lation are in progress and it is also envisaged to con- 
struct a volume-limited sample, spanning a reasonably 
large range in luminosity and mass, to test how much 
the Lx-Mtot relation given here is affected by being es- 
timated partly from a flux-limited sample. 

5.3. Mass Function 

5.3.1. Comparison to previous estimates 

Bahcall & Cen (1993) give a mass function constructed 

a) from optically selected clusters with masses deter- 
mined from the galaxy richness and b) from the clus- 
ter X-ray temperature function given by Henry & Ar- 
naud (1991). Very good agreement is found for masses 
determined within ta between the Bahcall & Cen and 
HIFLUGCS mass function (Fig. 13). 

White et al. (1993a) constrain the cluster abmidance 
by using published values for the abundance and median 
velocity dispersion of richness class R > 1 Abell clus- 
ters. It is not surprising that their density is significantly 
higher than the HIFLUGCS density since they have in- 
tentionally used conservative mass estimates, which are 
overestimates of the true cluster masses. 

Biviano et al. (1993) and Girardi et al. (1998) have 
determined the cluster mass function using optically se- 
lected cluster samples with masses determined from pub- 
lished line-of-sight velocity dispersions of cluster galaxies. 
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Fig. 13. Cumulative gravitational mass functions for 
HIFLUGCS using three different outer radii. Vertical 
error bars give the Poissonian errors. Horizontal bars 
indicate the individual bin sizes. Each bin contains 10 
clusters, apart from the highest mass bin, which contains 
13 clusters. The abundances from previous works are 
determined for Ma- 



At the completeness limit given by Girardi et al. (1998, 
triangle in Fig. 13) the cluster density given by Biviano 
et al. (1993) is about a factor of two higher than the 
density given by Girardi et al. (1998). The latter au- 
thors explain this by their on average 40 % smaller mass 
estimates due to an improved technique for removing in- 
terlopers and the use of a surface-correction term in the 
virial theorem. The value given by Girardi et al. itself 
lies significantly higher than the comparable HIFLUGCS 
density. The reason could lie in the fact that their opti- 
cally estimated masses arc in general slightly larger than 
the X-ray masses or that the external normalization for 
R > 1 {Ngx. > 50) clusters which they applied to their 
mass function is intrinsically higher than the normaliza- 
tion obtained through HIFLUGCS directly, or both. By 
comparing the mass estimates for a common subsample 
of 42 clusters it has been found that the virial masses 
determined by Girardi et al. (1998) are on average 25% 
larger than the masses determined in this work. This 
difference might be smaller if masses out to the Abell 
radius were compared, since common radii would have 
been used in this case. Increasing the X-ray masses ar- 
tificially by 25%, the diamonds in Fig. 13 shift towards 
higher masses, but the shift is too small to account for the 
difference to the triangle. The large scatter in the Ng^- 
Mtot diagram (Fig. 12) makes a reliable estimate of a best 
fit relation between these two quantities very difficult. 
Nevertheless, in order to get a rough idea of the mass 
within TA that corresponds to Ng^ = 50, we have per- 
formed fits using the minimization methods outlined ear- 
lier and find 5.1 < Ma{N^ = 50) < 8.8 x lO" h^^ Mq. 
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Note that this range is in agreement with the ranges 5- 
8 X 10^^ Mq and 5- 7 x lO" h'^ Mq given by Bah- 
call & Cen (1993) and Girardi et al. (1998), respectively, 
for A^gx = 50. This mass range corresponds to a cumu- 
lative number density obtained through HIFLUGCS in 
the range 1.7 < n(> Ma) < 8.7 x 10'^ hloMpc-^ . The 
external density estimate applied to normalize the Gi- 
rardi et al. mass function therefore is a factor 1.2-6.2 
higher than the corresponding estimate obtained here. 
It is therefore concluded that both effects (masses and 
normalization) are important but the latter factor is re- 
sponsible for a larger fraction of the discrepancy. Assum- 
ing both normalizations to be determined from samples 
that are highly complete and representative of the local 
universe this may indicate that either X-ray and opti- 
cal chistcrs arc drawn from different populations or that 
projection effects (e.g., line of sight galaxy overdensities, 
which do not form a bound structure in three dimen- 
sions) possibly bias optically determined normalizations 
high. 

Girardi & Giuricin (2000) recently extended the mass 
function to loose galaxy groups, finding n(> l.Sh'^Q x 
10" Mq) = 1.6 2.4 X IQ-^ hl^Mpc-^, which is outside 
the mass range wc can test. They find that the group 
mass function can be described by a smooth extension of 
the cluster mass function by Girardi et al. (1998). Con- 
sistently this abundance is higher than the abundance 
given by Bahcall & Cen (1993) at that mass scale. 

Carlberg et al. (1997) have compiled and partially 
rocstimatcd abundances for local cluster samples (Henry 
& Arnaud 1991; Mazure et al. 1996; White et al. 1993a; 
Eke et al. 1996) for comparison with higher redshift sam- 
ples (the ' X ' shows the density for a sample with higher 
mean redshift and therefore it cannot be compared di- 
rectly). Note that Borgani et al. (1999a) find that the 
reestimate of the mass limit of the Mazure et al. sample 
by Carlberg et al. (the '+' at n w 3 x IQ-^hloMpc-^) 
may lead to an underestimated mass limit. In general 
the comparison to the HIFLUGCS mass function shows 
better agreement than the abundances given by Girardi 
et al. (1998) and White et al. (1993a). 

The obvious importance of the definition of the clus- 
ter outer radius for the cluster mass function can be 
directly appreciated by noting the large differences be- 
tween the mass functions determined for HIFLUGCS for 
M500, M200, and Ma in Fig. 13, especially towards lower 
masses. Since for self similar clusters the mass scales with 
the third power of the characteristic radius (Sect. 5.2), 
determination of the mass within a characteristic over- 
density is the natural choice. We mainly give the formally 
determined Ma mass function for the comparison with 
previous mass functions and recall again that especially 
for the low mass systems the assumption of hydrostatic 
equilibrium may not be justified out to ta, and therefore 
our mass estimates of M500 and M200 should be consid- 
ered much more precise than the estimates for Ma. 



5.3.2. Comparison to predicted mass functions 

We use the standard formalism based on the Press- 
Schechter (PS) prescription to predict cluster mass func- 
tions for given cosmological models (see, e.g., Borgani 
et al. 1999b). To allow easier comparison with the 

theoretical literature on this subject in this paragraph 
hiao = /i5o/2 is used. The mass function is then given by 



dn{M) 
dM 



2 Po 5c{z) 
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exp 



'2a{Mf 



(6) 



(Press & Schechter 1974, Bond et al. 1991; see, e.g., 
Schuecker et al. 2001a for a compilation of published 
extensions of the PS mass function). Here M repre- 
sents the halo (cluster) virial mass and po = 2.7755 x 
10^^ rim /ifoo -^0 Mpc~^ is the present day mean mat- 
ter density. The linear overdensity computed at present 
6c{z) = 5l{z) D{0) D{z)^^ , where the linear overden- 
sity at the time of virialization, S^{z), is computed 
using the spherical collapse model summarized in Ki- 
tayama & Suto (1996), for = 1 using (A2) and for 
rim < 1 A r^k = using (A6,7); the linear growth factor 
dIz) =2.5n^E{z) J^{l + z')E{z')-^dz' and E{z) has 
been defined in Sect. 3.2. As mentioned earlier due to the 
low redshift range spanned by HIFLUGCS, the effect of 
a redshift correction is very small and we therefore set 
2: = for all calculations, unless noted otherwise. The 
variance of the cosmic mass density fluctuations 



o-(M)2 



, J^e+''T{kf\W{kR{M))\^dk 
' A;2+" T(fc)2 \W{k 8 /i^oo Mpc)|2 dk 



(7) 



where (Tg represents the amplitude of density fluctuations 
in spheres of radius 8 /i^qo Mpc. Recent measurements of 
the cosmic microwave background (CMB) anisotropies 
indicate that the primordial power spectral index, n, has 
a value close to 1 (e.g., Balbi et al. 2000; Jaffe et al. 
2001; Pryke et al. 2001; Wang et al. 2001; de Bernardis 
et al. 2001) and is therefore set to 1, unless noted other- 
wise. For the transfer function we use the fitting formula 
for Cold Dark Matter (CDM) cosmologies provided by 
Bardeen et al. (1986) for q{k) = fc/(r /iioo Mpc'^) 



[1- 



T{k) = T{q{k)) = ln(l - 
3.89g-h (16.1g)2 -t- (5.46g)3 



2.34^)7(2. 34g) 
-{6.71qrr'^\ (8) 



where the shape parameter is given by (modified to ac- 
count for a small normalized baryon density fib > 0, 
Sugiyama 1995) 



f2m h 



m 'noo 



2.7K 



exp 



1100 



0.5 On 



(9) 

The temperature of the CMB To = 2.726 K (Mather et al. 
1994) and Ob/i^oo = 0.0193 (Buries & Tytler 1998), 
for the latter equation and (9) /iioo = 0.71 has been 
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used (Mould et al. 2000). The comoving filter radius 

R{M) = [3M/(47rpo)]^/^ for the top hat filter function 
W{x) = 3 (sin a; — xcosx)/x^, which is adopted in this 
analysis, because the HIFLUGCS masses have been de- 
termined with a top hat filter, too.^ Since the PS recipe 
as outlined above assumes virial masses based on the 
spherical collapse model we use M200 as approximation 
to the virial mass (Sect. 3.2). 

Similarly to the work of Ikebe et al. (2001) the statisti- 
cal uncertainty in the mass determination is incorporated 
in the model mass function as 

dnjM) ^ 1 r dn{M') 
X exp / -aogM;-logM)A ^^^^ 

where ^Mtotiiog = 0.12 represents the logarithmic mean 
mass measurement uncertainty. Note that since HI- 
FLUGCS is flux- limited and not volume- limited the 
weighting has to be performed on the mass distribution, 
N(M)/dM, rather than on the mass function itself. The 
effect of this weighting on the model mass function is a 
slight amplitude increase at the high mass end. 

For the modeling to be independent of the precise 
knowledge of the Lx^-^^tot relation the quantitative com- 
parison has been performed using a standard proce- 
dure on the differential binned mass function given in 
Fig. 14 (rather than using a maximum likelihood ap- 
proach on the mass distribution). After identifying the 
crude position where is minimal in a large 

m cTs pa- 
rameter space region the values have been calculated 
in a fine grid of 200 by 200 Om-crg values in the range 
0.05 < < 0.26 and 0.65 < (7& < 1.30. A flat cos- 
mic geometry has been assumed, i.e. + f^A = 1- The 
cosmological constant enters the calculation only through 
(5c, however, and therefore has a negligible influence here. 
The minimum and statistical error ellipses for some stan- 
dard confidence levels (c.l.) are given in Fig. 15. The tight 
constraints obtained show that with HIFLUGCS we can 
go beyond determining an Vl^-a% relation and put limits 
on Om and erg individually. It is found that 

Q^ = 0.12to°4 and ag = 0.96t°;l^ (11) 

(90% c.l. statistical uncertainty for two interesting pa- 
rameters) , indic;ating a relatively low value for the density 
parameter. The large covered mass range, the specific re- 
gion in Om-cg parameter space, and the assumption of 
CDM cosmological models with given Hq, fib, and To 
allow to derive these tight constraints from a local cluster 



This approach follows the custom of disregarding the inconsistency 
of ijsing top hat masses while the PS mass function with the correct 
normalization has been derived for the sharp fc-space filter (Bond 
et al. 1991); see Schuecker et al. (2001a) for a generalization to 
more realistic filter functions. 



sample. For comparison for a given a% value the Om value 
which minimizes x^ is calculated. In the interval shown 
in Fig. 15 these pairs can then roughly be described by 
a straight line in log space given by 

c78 = 0.43n-°-38. (12) 

In Fig. 14 we also plot the best fit model mass func- 
tions for given fini = 0.5 and flm = 1.0 and one notes 
immediately that these value pairs give a poorer descrip- 
tion of the shape of the mass function. The slight dif- 
ferences between the best fit values for (Tg and the ones 
expected from (12) are caused by the fact that the sim- 
ple power law is only an approximation to the real shape 
of the error ellipse and weakly by the extrapolation in- 
volved. Previous estimates have generally yielded a com- 
bination of slightly higher Qm-CTg values (e.g., Peacock 
1999, Sect. 17.2). It has to be noted, however, that for 
instance in the important work of White et al. (1993a), 
who find (Tg w 0.57 the authors explicitly state 

that their estimates of ag are probably biased high due 
to their conservative mass estimates. 
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Fig. 14. — HIFLUGCS mass function compared to the 
best fit model mass function with Q,^ = 0.12 and erg = 
0.96 (solid line). Also shown are the best fit model mass 
fimctions for fixed ri,„ = 0.5 (^ (Jg = 0-60. dashed line) 
and fim = 1.0 (^ cTg = 0.46, dotted line). 

Before exploring sources of systematic errors let us 
evaluate our treatment of the shape parameter F. By 
applying the Bardeen et al. (1986) fitting formula to cal- 
culate the transfer function we have adopted the assump- 
tion of a CDM universe, i.e. Om — ^^b = ^cdm- Within 
this model V and Om are not independent. For a small 
r2b/^^m ratio F is given by (9). However, to test its in- 
fluence we have set F to a fixed value of 0.19 as im- 
plied by Large Scale Structure (LSS) measurements of 
galaxies and clusters (e.g., Szalay et al. 2001; Schuecker 
et al. 2001b, see also the discussion at the end of this 
Section). The best fit values for fim and erg change sig- 
nificantly, the Om value becoming smaller and the ug 
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value larger. The ilm-crs relation changes slightly to 
about as = OAOn^"-*^ for the range [0.04-0.30] 
(note that for this range relation (12) translates to about 
as = 0.39 0"" *^, the steepening being mainly due to a 
strong increase of as for the very lowest fim values) . The 
error ellipse in general becomes much longer when F is 
set to a fixed value, preventing the possibility to obtain 
tight individual constraints on and as, even though 
r^iji — 1.0 is still excluded at least at the 90% c.l. for 
all r > 0.03. The decreased sensitivity is due to the 
loss of the amplification of a change of flm through F (to 
first order oc F, eq. 9), since both a larger (smaller) 
r2,n and a larger (smaller) F result in steeper (less steep) 
model mass functions. Having found a dependence of 
the results on F we tried to constrain F directly from the 
observations, as done by previous authors (e.g., Borgani 
et al. 1999b). We have set F to 31 different fixed values 
in the range 0.01 < F < 0.50 and calculated a Xmin fo'" 
each of those by varying fi,„ and as ■ The result is shown 
in Fig. 16. The conclusion is that the data do not contain 
enough information to provide an independent constraint 
on F. Therefore we have taken advantage of the depen- 
dence of F on ilm within the CDM framework and have 
continued to calculate F with cq. 9. Notice that this pro- 
cedure requires additional knowledge about Hq and Qb- 
The influence of the uncertainty of these two parameters 
is tested along with various other possible systematic un- 
certainties in the following. It is important to note that 
the best fit value for Om is not determined solely through 
(9) but also directly through po in (6) and through the 
filter radius R{M). 

The quoted error ranges have been calculated from 
the procedure. Most of the remaining part of this 
Section is devoted to an identification and quantification 
of possible systematic uncertainties. The quantitative 
results are summarized in Tab. 5. 

Due to the large given ranges of several orders of mag- 
nitude in mass and especially density the values have 
been determined from comparison model/data naturally 
in log space. However, we have verified that the same cal- 
culation in linear space yields best fit values lying within 
the 68 % error ellipse. 

In Sect. 4.2 arguments have been given why we have 
used Vinax(-^x) for the determination of the mass func- 
tion. Nevertheless if V^ax(-^tot) is used instead (see 
Fig. 7), consistent results are obtained. Since in this 
case we want to estimate L from M the relation {L \ M) 
is the appropriate one (e.g., Isobe et al. 1990). Perform- 
ing a fit to a mass function constructed with 14nax(Aftot) 
results in best fit values fini = 0.14 and as = 0.86, which 
is consistent with the error range given in (11). 

The data point in Fig. 14 that may be affected most 
by cosmic variance is the one at the lowest mass, since 
the maximum search volume is smallest for the clusters 
in this bin. We therefore tested the sensitivity of the 



best fit results on this last point by ignoring it. The 
decrease in the covered mass range of course increases 
the resulting error ellipse, but the best fit values vary 
only within the (smaller) 68% error ellipse. It may be 
worth noting that leaving out the highest mass bin or 
leaving out the highest and lowest mass bin changes the 
best fit values only within the 90 % error ellipse. 

Simulations seem to indicate that the mass den- 
sity profile of virialized halos follows the NFW profile 
(Navarro et al. 1996; Navarro et al. 1997, but see, e.g., 
Kravtsov et al. 1998). As test therefore gravitational 
masses have been recalculated assuming the dark matter 
density to follow an NFW profile based on the measured 
values for f3 and Vc- The approximation Vg = rc/0.22 
and b = (3/0.9 given by Makino et al. (1998) has been 
used combined with their eqs. 7 and 9. r2oo and M200 
have been redetermined based on this model and the re- 
sulting mass function has been compared to model mass 
functions. On average it is found that the NFW masses 
are a factor of 0.90 lower, the difference being smaller for 
high mass clusters than for low mass ones. The values 
for ri,ji and erg for which is minimal lie within the 
90% statistical error ellipse shown in Fig. 15. Therefore 
the choice of mass calculation does not affect the results 
significantly. 

As shown in Sect. 3.2 it is possible that the assumption 
of isothermality leads to an overestimate of the cluster 
masses on average. Therefore the robustness of the re- 
sults has been tested by irmltiplying the isothermal clus- 
ter masses by 0.80 and recalculating the minimum. As 
expected values for both and as are found to be lower. 
But the new minimum is contained well within the error 
ranges given in (11). On the other hand in Sect. 5.3.1 it 
has been shown by comparison with optical virial mass 
measurements that masses could possibly be underesti- 
mated. After increasing all cluster masses by 25 % a fit 
shows that both and as become slightly larger, but 
again they lie well within the range (11). Therefore these 
tests indicate that the constraints obtained here are fairly 
insensitive against systematic uncertainties in the mass 
measurements. Note that systematic effects that show 
a mass dependence have a higher potential to affect the 
individual best fit values of Om and erg than effects which 
simply shift all data points by about the same amount in 
the same direction. However, no indications for a mass 
dependence of the differences between the above mass 
estimations have been found. 

It has been mentioned that masses calculated within 
rsoo are less affected by systematic uncertainties because 
less extrapolation beyond rx is needed and because ac- 
cording to simulations the assumption of hydrostatic 
equilibrium can be more safely applied. Nevertheless we 
have used M200 for the comparison to predicted mass 
functions because it is a better approximation to the 
cluster virial mass. However, for comparison the best fit 
values have also been determined using M500, knowing 
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Fig. 15. Statistical confidence contours for the pro- 
cedure. The cross indicates the position of the mini- 
mum, Xmin- ElUpses indicate the 68%, 90%, 95%, and 
99 % confidence levels for two interesting parameters, i.e. 
Ax' = X' - Xmin = 2.30, 4.61, 6.17, and 9.21, respec- 
tively. 

that virial masses arc probably underestimated this way. 
The result is that the best fit value for fim changes only 
marginally. The best fit (jg value is slightly lower than 
allowed by the statistical error. On the other hand the 
spherical collapse model implies that one may need to ex- 
trapolate even further than r2oo for low density universes. 
For instance in a flat universe Om =0.2 implies that the 
virial radius is close to rs?. Therefore a mass function 
has been constructed using Mgy as cluster masses, know- 
ing that these masses are likely to be rather uncertain. 
Nevertheless the resulting best fit flm value again varies 
only insignificantly whereas the erg value becomes signif- 
icantly larger. Note that for both tests the best fit value 
for Om hardly changes. 

As mentioned in Sect. 3.2 it has been shown that the 
small range of low redshifts covered here ensures that 
no redshift corrections need to be applied. Nevertheless 
we have tested whether or not the best fit parameter 
values change if M20Q is calculated using pc = Pc{z) for 
the extreme (strong evolution) case (flm = l;f^A = 0), 
i.e., Pc = 4.6975 x 10~^° (z + l)^gcm~^ for each cluster 
redshift. The model mass function is then calculated for 
the mean redshift of HIFLUGCS, {z) = 0.05, using the 
formulae outlined in the beginning of this Section. We 
have found that within our grid the best fit values do 
not change at all and also the error ellipses arc almost 
not affected, thereby confirming that the application of 
redshift corrections does not affect the results. 

The value Hq = 71kms^iMpc"i has been adopted 
for the calculation of model mass functions based on 
the recent combination of constraints obtained using 
the Hubble Space Telescope (Mould et al. 2000). Set- 
ting the Hubble constant to their lower limit, Hq = 
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Fig. 16. Minimal x' values calculated by varying f7m 
and as for a variety of fixed values of F instead of using 
eq. 9. The dashed line represents the 68% c.l. statisti- 
cal uncertainty for one interesting parameter. This line 
actually underestimates the uncertainty since the proce- 
dure involves three interesting parameters. Obviously F 
cannot be constrained significantly. 

65 kms~^ Mpc~^, does not affect the best fit parameter 
values significantly. Using even Hq — 60kms~^Mpc~^ 
changes the results for fim and erg only well within the 
68 % error ellipse. In Tab. 5 results have also been listed 
for Hq = 50 and 80kms~^Mpc^^ showing that even 
variations within this range change the fim and erg values 
only within their statistical uncertainties given in (11). 
Therefore our constraints on these cosmological parame- 
ters do not depend significantly on the specific choice of 
Ho. 

A baryon density based on measurements of the pri- 
mordial deuterium to hydrogen ratio combined with stan- 
dard nucleosynthesis theory has been used in this work. 
The used value also agrees with recent CMB measure- 
ments (e.g., Pryke et al. 2001). However, we have also 
calculated the minimum for the values fib h^oo = 0.0100 
and 0.0300 finding that the position of the minimum in 
each case does not deviate significantly from our original 
result. 

The primordial power spectral index has been set to 
n = 1 based on measurements of the CMB fluctuations. 
The influence of a change in this parameter has been esti- 
mated by calculating the minimum for the values n = 0.8 
and n = 1.2. The differences in the best fit parame- 
ters are smaller than the statistical uncertainties given 
in (11). 

Since we have found that for the estimate of the sta- 
tistical errors we need to explore ranges il,n < 0.1 for 
fib ~ 0.04 we regarded it necessary to check if the approx- 
imation to the transfer function as given in the beginning 
of this Section is still applicable. Recently Eiscnstcin & 
Hu (1998) derived a fitting function, that includes, e.g.. 
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also the oscillations induced by the baryons, which gives 
a better description of transfer functions computed with 
CMBFAST (Seljak & Zaldarriaga 1996) than fitting func- 
tions for zero or small baryon contribution to the total 
matter density derived previously. Therefore we incorpo- 
rated this improved version of an analytic transfer func- 
tion in the procedure. We have found that within 
our grid there is only a very weak shift of the minimum, 
the choice of the Bardeen et al. (1986) fitting function 
combined with the shape parameter given by Sugiyama 
(1995) therefore seems to be accurate enough for our pur- 
poses. However, the confidence contours towards low flm 
are getting compressed when the Eisenstein & Hu (1998) 
fitting function is used, thereby slightly decreasing the 
area of the error ellipse for a given confidence level. Since 
we regard the latter statistical error ellipse as more real- 
istic, we show this one in Fig. 15. 

We have also tested whether or not the recently found 
deviations of the PS formalism compared to large A^- 
body simulations (e.g., Governato et al. 1999; Jenkins 
et al. 2001) have a significant influence on the results ob- 
tained here. We have compared the best fit PS model 
(Qm = 0.12, (7s = 0.96) to the model obtained using 
the 'universal' mass function (fit to A'^-body simulations, 
Jenkins et al. 2001) for the same parameter values. These 
two models agree well for M < lO^^ h^^ Mq. The differ- 
ences become larger than the size of the Poissonian error 
bars (Fig. 14) for M > 2 x IO^^so^Mq, in the sense 
that the Jenkins et al. mass function predicts higher 
cluster abundances than PS. For larger values of ri,„ the 
differences become comparable to the size of the error 
bars at lower masses, e.g., for fl^ = erg = 0.5 around 
M ^ 5 X 10^** h'^Q Mq. To estimate the influence of these 
differences on the best fit values derived using the PS 
mass function, we adjusted the parameter values of the 
Jenkins et al. model to reproduce the PS mass func- 
tion, finding — 0.15 and a"8 = 0.86. The value for 
fim becomes slightly larger but the combination of both 
values is still contained within the 90 % error ellipse. We 
therefore conclude that the differences between the model 
mass functions do not significantly affect the interpreta- 
tion of the HIFLUGCS mass function. Moreover we re- 
gard this test as confirmation of the validity of the PS 
mass function for the accuracy needed here. 

Almost all identified systematics give rise to smaller 
uncertainties than the statistical uncertainties (eq. 11). 
However, it is not impossible that several systematics 
combined have a significant effect. In order to determine 
the highest fim value allowed by a conspiracy of all iden- 
tified uncertainties, all positive Afim values from Tab. 5 
have been added to the positive statistical uncertainty for 
rim and this sum has been added to the l^m value which 
minimizes (^m = 0.12). In this worst case procedure 
the upper limit 

Clm < 0.31 (13) 



has been found. 

Using = 0.19 ± 0.08 h'^Q^'^ ^m, as derived from the 
mean gas mass fraction within r2oo for the 106 clusters 
in the enlarged sample, results in model mass functions 
very similar to the ones calculated using the baryon frac- 
tion given by Buries & Tytler (1998). It is therefore not 
surprising that the best fit values for fim and ag vary 
only well within the 68 % error ellipse if the former 
determination is used. Moreover it is worth noting that 
combining these two measurements of fib gives further 
evidence for a low value for fim by yielding an estimate 
= OMtoAO using Hq = Tlkms^^Mpc-^ (Mould 
et al. 2000), where the error has been determined from 
the standard deviation of fib given above. This value for 
Sim is an upper limit since baryons, e.g., contained in the 
cluster galaxies or forming part of the dark matter have 
been neglected. A low value for fim has been indicated by 
this method for smaller cluster samples by several works 
previously (e.g.. White & Frenk 1991; Briel et al. 1992; 
Bohringer 1993; White et al. 1993b; White & Fabian 
1995; Ettori & Fabian 1999; but see Sadat & Blanchard 
2001). One has to keep in mind, however, that this es- 
timate extrapolates the gas fraction from cluster scale 
to cosmic scales. For the clusters in our sample we have 
found that the gas fraction is not constant but varies with 
radius and cluster mass (Reiprich & Bohringer 1999b), 
therefore further observational tests of this assumption 
may be useful. 

Let us briefly discuss the role of F again. The con- 
straints (11) together with hwo = 0.71, fly, /ifoo = 0.0193, 
and eq. 9 imply F = OMtoit 

which is significantly 
lower than values for T determined through LSS stud- 
ies of galaxies (e.g., Szalay et al. 2001 have found F = 
0.188 ± 0.04) or clusters (e.g., Schuecker et al. 2001b 
have found F = 0.195 ± 0.055). Adopting the worst case 
upper limit on ri,„ found above yields an upper limit 
F < 0.18 consistent with the LSS results. However, 
also many previous studies of cluster abundances have 
yielded low values for F. For instance Eke et al. (1998) 
have fomid F w 0.09 ± 0.08; Henry (2000) has found 
F « 0.05 ± 0.22; Oukbir & Arnaud (2001) have found 
F = 0.05fi-°-^i+°-^<5"» which, e.g., for fl^ = 0.3 yields 
F w 0.09; judging from their Fig. 3 and their best fit val- 
ues for a spatially flat universe Borgani et al. (2001) have 
found F « 0.1, however, not significantly constrained. It 
is unclear whether these values indicate a significant dis- 
crepancy but it is interesting to note that cluster abun- 
dances seem to favor lower values for F than LSS studies 
do. If significant this could possibly imply either uniden- 
tified systematic errors in one or both measurements - 
from the cluster abundance side one could for instance 
imagine a higher merger rate for more massive clusters 
resulting in overestimated X-ray luminosities, tempera- 
tures, and masses for the massive clusters resulting in 
less steep Lx, Tx, and Mtot functions and therefore in 
a decreased value for F (and possibly fim) - or that the 
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Table 5 
Systematic uncertainties 





0-8 




Ao-8 


Test 


0.116 


0.928 


±0.000 


-0.026 


(1) 


0.137 


0.865 


+0.021 


-0.094 


(2) 


0.150 


0.897 


+0.034 


-0.062 


(3) 


0.092 


1.082 


-0.024 


+0.124 


(4) 


0.083 


1.050 


-0.034 


+0.091 


(5) 


0.150 


0.860 


+0.034 


-0.099 


(6) 


0.107 


0.907 


-0.009 


-0.052 


(7) 


0.132 


1.004 


+0.016 


+0.046 


(8) 


0.106 


0.826 


-0.011 


-0.133 


(9) 


0.108 


1.242 


-0.008 


+0.281 


(10) 


0.116 


0.959 


±0.000 


±0.000 


(11) 


0.118 


0.956 


+0.002 


-0.003 


(12) 


0.153 


0.861 


+0.037 


-0.098 


(13) 


0.108 


0.991 


-0.008 


+0.033 


(14) 


0.106 


1.001 


-0.011 


+0.042 


(15) 


0.128 


0.923 


+0.012 


-0.036 


(16) 


0.148 


0.874 


+0.032 


-0.085 


(17) 


0.094 


1.056 


-0.022 


+0.098 


(18) 



Note. — Sec text for a more detailed de- 
scription of the tests. The differences to the 
best fit results have been been calculated with 
more digits than given here. 

References. — (1) calculated in lin- 
ear space. (2) Vmax(-Mtot) used with rela- 
tion (L I M) for HIFLUGCS (Tab. 8). (3) 
Lowest mass bin ignored. (4) Highest mass 
bin ignored. (5) i"2ao and M200 calculated 
assuming NFW profile. (6) Universal mass 
function. (7) M2m) reduced by factor 0.80. 
(8) M200 increased by factor 1.25. (9) Afsoo 
used. (10) MsT used. (11) Redshift correction 
applied for cluster mass determination and 
model mass function calculation. (12) Trans- 
fer function calculated using Eisenstein & Hu 
fitting formula. (13) hum set to 0.5. (14) 
/iioo set to 0.8. (15) set to 0.0100 ^^qq. 
(16) Ob set to 0.0300 (17) n set to 0.8. 
(18) n set to 1.2. 
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Table 6 

Examples of previous results for given fl, 



as Relation Ref. 



0.30 


1.12 


<^8 


= 0.570^°-^^ 


1 


0.30 


1.16 


0"8 


^ ^„ „-0.59+0.16f2m-0.06n?, 
= 0.60 Sim 


2 


0.30 


0.93 


0-8 


= 0.52a-''-52+o.i3n„ 


3 


0.30 


1.01 


<^8 


= o.6on^i°'*«+°°8«n> 


4 


0.30 


0.78 






5 


0.30 


0.99 


<^8 


= 0.56 0-°*'' 


6 


0.30 


0.96 




= o.58n-o''^+°-i«"° 


7 


0.30 


0.96 






8 


0.30 


1.02 


0"8 


= 0.495 


9 


0.30 


0.91 


0-8 


r, ,.„r^-0.57+1.45nm-3.48n?,+3.77n5,-1.49ni 

= 0.59S2ni m m m 


10 


0.30 


0.68'' 


as 


= 0.43 0^°-^* 


11 



Tor consistency the iT8 value has been calculated using the relation; the 
directly determined best fit value for Hm = 0.30 is given by its = 0.72. 



References.— (1) White et al. 1993a. (2) Viana & Liddle 1996. (3) 
Eke et al. 1996. (4) Girardi et al. 1998. (5) Markevitch 1998. (6) Viana 
& Liddle 1999. (7) Borgani et al. 1999b. (8) Blanchard et al. 2000. (9) 
Pierpaoli et al. 2001. (10) Oukbir & Arnaud 2001. (11) This work. 



description of the power spectrum using a single T pa- 
rameter on all scales is not accurate (e.g., Eisenstein & 
Hu 1998; Saslaw 2000, Sect. 33.3.3). 

Even though we made a conservative estimate by ne- 
glecting the possible presence of gas temperature gra- 
dients, previous estimates obtained from cluster abun- 
dances generally yielded higher values for Om and as 
(e.g., White et al. 1993a; Girardi et al. 1998). For the 
latter two works one could have expected this already 
from Fig. 13 and possible reasons have been discussed 
in Sect. 5.3.1. Recently a number of constraints on 
and (Tg have been obtained using different kinds of dis- 
tribution functions assuming relations between the mea- 
sured quantities and mass. Some of them incorporated 
additional dynamical information from the evolution of 
the corresponding distribution function. Especially intr- 
acluster gas temperature functions have been constructed 
and used frequently. In Tab. 6 some results, which al- 
lowed a simple comparison for a given = 0.3, are 
listed. If there was a choice, results assuming flat cosmo- 
logical models and isothermal clusters are quoted. This 
Table is not a complete summary of recent results but 
merely shows that previous results yielded higher values 
for (Ts than obtained here from a newly constructed sam- 
ple (even using Mgr instead of M200 would increase the 
as value only up to 0.82). The reasons may be manifold. 
Some have been discussed above. Other reasons may in- 
clude different sources of temperature estimates. For in- 
stance Pierpaoli et al. (2001) used cooling flow corrected 



temperatures given by White (2000), which are generally 

higher than temperatures derived previously. Using the- 
oretical/simulated Mtot-T'gas relations (e.g., Evrard et al. 
1996; Bryan & Norman 1998) to connect observed tem- 
perature functions with theoretical mass functions will 
on average also result in higher values for as , since these 
relations have been shown to have higher normalizations 
than observed X-ray mass-temperature relations (e.g., 
Horner et al. 1999; Nevalainen et al. 2000; Finoguenov 
et al. 2001). Assuming low temperature groups to be 
affected by non gravitational heat input, temperature 
functions, unlike X-ray mass functions which incorpo- 
rate the additional information of the gas density profile, 
would show an additional artificial increase towards low 
temperatures. If not taken into account, e.g., by use 
of a modified Mtot^T'gas relation, this effect would re- 
sult in artificially increased values for 17m (see Fig. 14). 
Also the completeness and representativeness of the used 
cluster samples obviously play an important role for the 
^^m-cTs determination. Here advantage has been taken 
of the currently best available local sample in the sense 
of homogeneous X-ray selection, completeness, size, and 
availability of high quality observations. The results ob- 
tained here are, however, in good agreement with the 
results from the power spectral analysis of the 452 RE- 
FLEX clusters. Schuecker et al. (2001b) find for a given 
ACDM model (f7i„ = 0.3) that as = 0.7 represents the 
data well, which is very close to as = 0.68 expected us- 
ing relation (12) found here. Moreover the (l-cr) range 
0.17 < fim < 0.37 (using h = 0.71 in their eq. 18) quoted 
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for Qm directly is also consistent with the 90 % range de- 
termined hero. Furthermore Ikobc ct al. (2001), who an- 
alyzed the HIFLUGCS temperature function using tem- 
peratures from homogeneously reanalyzed ASCA data, 
find an = O.lStJ] ;?^ and = OMt'^H (90% c.l. statis- 
tical uncertainty; assuming an open cosmology) by com- 
parison with Press-Schechter models, which is in good 
agreement with our results. Moreover, after this paper 
had been submitted another paper appeared on astro-ph 
(Borgani ct al. 2001), where the authors also find Om-crg 
values that agree with values calculated using relation 
(12) presented here. 

5.3.3. Mass function estimateA using Ly^ Aftot relation 

To show consistency and the power of the Lx-Mtot 
relation, we have also performed fits to 'mass' functions, 
where masses have been estimated from the measured 
X-ray luminosity. Relations for the flux-limitcd sample 
have been used to get the best mass estimate for the clus- 
ter luminosities included in HIFLUGCS. Mass functions 
for the two extreme relations {M \ L) with a = 1.538 and 
(L I M) with a — 1.310 are shown in Fig. 17. First of all 
one notes the fairly good agreement between the three 
mass functions. In detail the differences between the two 
mass functions estimated from different luminosity-mass 
relations can be understood by considering that at the 
low luminosity end the steeper relation predicts a higher 
mass for a given luminosity than the shallower relation, 
resulting in a shift towards higher masses of the mass 
function. At the high mass side the effect is opposite, 
resulting in a shift towards lower masses for the steeper 
relation. On average the points for the steeper relation lie 
higher which is caused by the fact that a steeper relation 
results in a smaller dM on average, which gives rise to an 
increased dn/dM. The differences to the mass function 
calculated using the measured masses are again under- 
stood by a similar comparison and are partially caused by 
a possible deviation of the shape of the Lx- Aftot relation 
from a pure power law. Despite these small differences 
performing an actual fit^ results in the (fimj cs) values 
(0.14, 0.85) for {L \ M) and (0.22, 0.74) for (M | L). The 
first case is consistent with the error range given in (11) 
and in the second case the 90 % statistical error ellipse 
overlaps with the 90 % ellipse in Fig. 15. From the above 
and Fig. 14 it is clear that using steeper luminosity-mass 
relations results in higher values for ri,„ and lower values 
for (J8. Here we want to estimate M from L and there- 
fore (M I L) is the appropriate relation to use. This test 
shows that with a comparatively easy to obtain X-ray 
luminosity function of a statistical cluster sample and 
with the knowledge of the empirical Lx--Mtot relation 

^For the fit the corresponding scatter in logM for the two relations 
(L I M) and (M | L), uiogMtot = 0-22 and 0.21, respectively 
(Tab. 12), has replaced the mass measurement error, CTMtot.log = 
0.12, in (10). 



(even if approximated as simple power law) and its scat- 
ter as presented here useful constraints on cosmological 
parameters can be set by construction of a 'quasi mass 
function'. 
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Fig. 17. HIFLUGCS mass function (filled circles) 
compared to 'mass' functions estimated using measured 
luminosities and luminosity-mass relations (open sym- 
bols). Squares have been calculated using the {M \ L) 
relation and triangles using the {L | M) relation for HI- 
FLUGCS clusters. 



5.4. Total Gravitating Mass in Clusters 

To estimate the fraction of the gravitational mass den- 
sity relative to the critical density contained in galaxy 

clusters, ^^Cluster(> Mtot,min) = l/Pc J^^^^ ^.^ Mtot 0(Mtot) dMtot, 

the individual cluster masses divided by the correspond- 
ing maximum search volumes have been summed up 

for HIFLUGCS, i.e., ^duster = l/Pc A^200,i/Vmax,i- 

Note that the determination of Ociuster is independent 
of the Hubble constant. The cumulative diagram for 
f^ciuster(> M200) is shown in Fig. 18. In order to per- 
form a conservative error estimate, HIFLUGCS has been 
split into two parts with bu > +20 deg and 611 < —20 deg, 
and the results for these subsamples are also shown in 
the Figure. This estimate is conservative because HI- 
FLUGCS is about twice as large as each subsample. 
Taking the second and third lowest mass clusters to- 
gether with the maximum mass range given by their 
individual uncertainties, we obtain 

^^Cluster = 0.012t°°°i (14) 

for masses larger than 6.4^Q g x 10^"^ h'^^ Mq, i.e., the 
total gravitating mass contained within the virial radius 
of clusters amounts only to 1.2^0 4 percent of the total 
mass in a critical density universe. Combined with our 
best estimate Om = 0.12 this implies that about 90% 
of the total mass in the universe resides outside virial- 
ized cluster regions above the given minimum mass. If 
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galaxies trace mass it also follows that by far most of 
the galaxies do not sit in clusters. This result is consis- 
tent with the general presumption that clusters are rare 
objects, rare peaks in the density distribution field. 




M2OO [hso"' Mq] 



Fig. 18. — Mass density contained in galaxy clusters 
as a function of minimum mass. Filled circles indicate 
the complete HIFLUGCS, open triangles indicate the 34 
clusters north of the galactic plane, and open diamonds 
the 29 clusters at southern galactic latitudes included in 
HIFLUGCS. 

Comparing the diagram to the mass fraction Ociuster = 
0.0281 

aoo8 clusters with masses larger than 2 x 
10" h^Q Mq given by Fukugita et al. (1998), based on the 
mass function determined by Bahcall & Cen (1993), one 
finds that their estimate is a factor 4-5 higher. However, 
the Bahcall & Cen mass function is given for Ma and we 
get a consistent result if we calculate cluster using our 
formally determined cluster masses within ta- It needs 
to be pointed out that at Ma ~ 2 x 10^^ h^^^ Mq we find 
that the typical virial radius is ~ 1 h^Q Mpc and a mass 
determination at 3 h'^Q Mpc based on the assumption of 
virial equilibrium may therefore be rather uncertain and 
possibly leading to overestimates of riciuster- This be- 
comes more crucial if mass functions for Ma are extrap- 
olated even down to galaxy masses. This way Fukugita 
et al. (1998) find Ocroup = 0.12 ± 0.02 within the mass 
range 2 x 10^^ 2 x 10^^ hT^Q Mq, which, compared to our 
results from the previous Section, would account already 
for almost all mass in the universe. 

6. Summary 

An X-ray selected and X-ray flux-limited sample com- 
prising the 63 X-ray brightest galaxy clusters in the sky 
(excluding the galactic band, called HIFLUGCS) has 
been constructed based on the ROSAT All-Sky Survey. 
The flux limit has been set at 2 x 10"^-^ ergss"^ cm~^ in 
the energy band 0.1 — 2.4 keV. It has been shown that 
a high completeness is indicated by several tests. Due 



to the high flux limit this sample can be used for a vari- 
ety of applications requiring a statistical cluster sample 
without any corrections to the effective survey volume. 

Mainly high quality pointed observations have been 
used to determine fluxes and physical cluster parame- 
ters. It has been shown that a tight correlation exists 
between the X-ray luminosity and the gravitational mass 
using HIFLUGCS and an extended sample of 106 galaxy 
clusters. The relation and its scatter have been quanti- 
fied using different fitting methods. A comparison to 
theoretical and numerical predictions shows an overall 
agreement. This relation may be directly applied in large 
X-ray cluster surveys or dark matter simulations for con- 
versions between X-ray luminosity and gravitating mass. 

Using HIFLUGCS the gravitational mass function 
has been determined for the mass interval 3.5 x 10^^ < 
M2Qa < 5.2 X 10^^ h'^Q Mq. Comparison with Press- 
Schechter mass functions assuming CDM power spectra 
has yielded tight constraints on cosmological parameters. 
The large covered mass range has allowed to put con- 
straints on relevant parameters individually. Speciflcally 
we have found = 0.12t°;|]^ and ag = OMto]_l (90% 
c.l. statistical uncertainty). Various tests for systematic 
uncertainties have been performed, including compari- 
son of the Press-Schechter mass function with the most 
recent results from large A^-body simulations, almost al- 
ways yielding uncertainties smaller than the statistical 
tmcertaintics. Combining all identified systematic uncer- 
tainties in a worst case scenario results in an upper limit 
fim < 0.31. For comparison we have also determined the 
best fit Om values for fixed erg values yielding the relation 
a8=0.43n-0-38. 

The mass function has been integrated to obtain the 
fraction of the total gravitating mass in the universe 
contained in galaxy clusters. Normalized to the critical 
density we have found r^ciustcr = 0.012^0.004 ^r cluster 
masses larger than 6.4^0.6 ^ 10^^ ^Q- With the value 
for Om determined here this implies that about 90 % of 
the mass in the universe resides outside virialized cluster 
regions. 
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Table 7 
Fit parameter values 



Fit 


a 


Aa 


A 


AA 


BCES(L 1 M) 


1.496 


0.089 


-17.741 


1.320 


bootstrap 


1.462 


0.089 


-17.238 


1.327 


BCES(M 1 L) 


1.652 


0.085 


-20.055 


1.261 


bootstrap 


1.672 


0.086 


-20.357 


1.278 


BCES-Bisector 


1.571 


0.083 


-18.857 


1.237 


bootstrap 


1.562 


0.083 


-18.717 


1.239 


BCES-Orthogonal 


1.606 


0.086 


-19.375 


1.283 


bootstrap 


1.609 


0.088 


-19.419 


1.308 



Note. — Best fit parameter values and standard devia- 
tions for the extended sample (106 clusters) for a fit of the 
form given in eq. 4. The rows denoted 'bootstrap' give the 
results obtained for 10 000 bootstrap resamplings. 



Table 8 
Fit parameter values 



Fit 


a 


Aa 


A 


AA 


BCES(L 1 M) 


1.310 


0.103 


-14.935 


1.526 


bootstrap 


1.256 


0.103 


-14.146 


1.531 


BCES(M 1 L) 


1.538 


0.105 


-18.320 


1.568 


bootstrap 


1.584 


0.113 


-18.995 


1.681 


BCES-Bisector 


1.418 


0.097 


-16.536 


1.434 


bootstrap 


1.407 


0.096 


-16.368 


1.427 


BCES-Orthogonal 


1.460 


0.105 


-17.157 


1.559 


bootstrap 


1.468 


0.110 


-17.274 


1.633 



Note. — Same as Tab. 7 but for the purely flux-limited 
sample (63 clusters). 



Table 9 
Fit parameter values 



Fit 


a 


Aa 


A 


AA 


BCES(L I M) 


1.756 


0.091 


-21.304 


1.350 


bootstrap 


1.719 


0.090 


-20.746 


1.338 


BCES(M I L) 


1.860 


0.084 


-22.836 


1.246 


bootstrap 


1.881 


0.084 


-23.144 


1.250 


BCES-Bisector 


1.807 


0.084 


-22.053 


1.251 


bootstrap 


1.797 


0.084 


-21.899 


1.241 


BCES-Orthogonal 


1.835 


0.085 


-22.473 


1.260 


bootstrap 


1.841 


0.085 


-22.563 


1.270 



Note. — Same as Tab. 7 but for LboI- 
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Table 10 
Fit parameter values 



Fit 


a 


Aa 


A 


AA 


BCES(L 1 M) 


1.504 


0.089 


-17.545 


1.298 


bootstrap 


1.469 


0.089 


-17.042 


1.300 


BCES{M 1 L) 


1.652 


0.086 


-19.708 


1.254 


bootstrap 


1.671 


0.086 


-19.992 


1.260 


BCES-Bisector 


1.575 


0.084 


-18.590 


1.228 


bootstrap 


1.565 


0.083 


-18.445 


1.224 


BCES-Orthogonal 


1.609 


0.087 


-19.075 


1.274 


bootstrap 


1.611 


0.088 


-19.109 


1.290 



Note. — Same as Tab. 7 but for M500. 



Table 11 
Fit parameter values 



Fit 


a 


Aa 


A 


AA 


BCES(L 1 M) 


2.683 


0.215 


-35.675 


3.219 


bootstrap 


2.606 


0.218 


-34.525 


3.265 


BCES(M 1 L) 


2.488 


0.127 


-32.761 


1.902 


bootstrap 


2.520 


0.130 


-33.239 


1.946 


BCES-Bisector 


2.583 


0.154 


-34.170 


2.297 


bootstrap 


2.560 


0.157 


-33.840 


2.347 


BCES-Orthogonal 


2.513 


0.128 


-33.137 


1.915 


bootstrap 


2.530 


0.132 


-33.388 


1.975 



Note. — Same as Tab. 7 but for Ma for comparison (the 
same relative mass errors as for M200 have been assumed 
here) . 



Table 12 
Measured scatter 



Scatter (L \ M) (M | L) Bisector 



106 clusters included in the extended sample. 



<^logMtot 0.21 0.21 0.21 

o-logix 0-31 0.34 0.32 

(^logL/M 0.17 0.18 0.17 



63 clusters included in HIFLUGCS. 



'^logMto, 0.22 0.21 0.21 

(JiogLx 0.29 0.32 0.30 



Note. — Scatter measured for different rela- 
tions. 
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